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A  general  method  is  presented  for  calculating  state-to-state 
cross  sections  for  photodissociation  of  polyatomic  molecules  by 
visible  and  UV  radiation.  The  method  also  provides  branching  ratios 
by  self-consistently  coupling  nuclear  motion  to  transitions  between 
several  electronic  surfaces.  Cross  sections  are  calculated  using  a 
transition  integral  between  the  initial  ground  state  wavefunction  of 
the  molecule  and  a  final  wavefunction  for  the  relevant  states  of  the 
product  fragments.  The  final  state  nuclear  motion  wavefunction  with 
incoming  wave  boundary  conditions  is  generated  in  an  eikonal  form 
along  classical  trajectories  which  follow  an  effective  potential  and 
contains  the  exponential  of  a  common  action  integral  and  a  pre- 
exponential  given  by  an  amplitude  matrix  and  a  Jacobian  describing  the 
divergence  of  the  trajectories.   In  a  new  development,  the  Jacobian  of 
the  variable  transformation  is  generated  exactly  along  a  classical 


VI 


trajectory  without  requiring  the  simultaneous  integration  of  adjacent 
trajectories.  Two  methods  are  presented  for  satisfying  the  asymptotic 
conditions  of  the  eikonal  wavef unction. 

In  model  calculations  on  CH-I  the  dynamics  are  described  by  two 
degrees  of  freedom — the  relative  position  of  I  to  CH-  and  the  umbrella 

mode  of  CH-.   The  first  calculation  involves  a  transition  to  a  single 

3 
dissociative  excited  Q^   potential  energy  surface.  The  total  and 

partial  cross  sections  for  an  initial  zero  vibrational  level  over  a 

range  of  photon  energies  are  found  to  be  in  excellent  agreement  with 

results  in  the  literature  from  an  exact  wavepacket  calculation  using 

the  same  empirical  potential  energy  surfaces.  Results  are  also 

presented  for  excitation  from  the  first  three  excited  vibrational 

levels  of  the  ground  electronic  state.  The  second  calculation 

involves  an  excitation  to  two  coupled  potential  surfaces  leading  to 

2         *  2 
CH-  and  I(  Pq/o)  or  I  (  ?..j).      Quantum  yields  of  .84  and  .69  are 

2 
calculated  for  I*(  P-i/o)  production  at  photon  wavelengths  of  266nm  and 

2 
248nm  respectively.  The  cross-section  curves  for  I(  Po/o)  production 

peak  at  higher  energies  than  the  results  from  a  coupled  channels 

calculation.   Possible  reasons  for  the  discrepancy  are  presented.   The 

method  is  also  applicable  to  large  polyatomic  molecules. 
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CHAPTER  I 
INTRODUCTION 


This  study  addresses  the  problem  of  formulating  and  testing  a 
computationally  efficient  and  accurate  theory  for  treating  the 
photodissociation  of  polyatomic  molecules  with  many  internal  degrees 
of  freedom  and  many  electronic  states.  The  photodissociation  event  of 
interest  will  be  limited  to  cases  where  incident  radiation  in  the 
ultraviolet  to  visible  range  is  absorbed  by  a  molecule  causing  a 
transition  to  a  dissociative  state  which  gives  a  single  atom  and  a 
molecular  fragment  as  products. 

There  are  already  many  theories  of  photodissociation.   They 
include  quantum  mechanical  approaches  [Heather  and  Light,  1983; 
Balint-Kurti  and  Shapiro, 1981;  Clary,  in  press]  as  well  as  classical 
and  semiclassical  trajectory  approaches  [Atabek  et  al.,  1976;  Mukamel 
and  Jortner,  1976;  Billing  and  Jolicard,  1984;  Gray  and  Child,  1984; 
Henriksen,  1985].   Lee  and  Heller  have  also  developed  a  theory  for  an 
exact  time  dependent  approach  using  wave  packets  [Lee  and  Heller, 
1982].   An  exact  approach  in  principle  gives  results  which  are  more 
reliable  than  a  semiclassical  method.   However,  for  large  molecules 
with  many  internal  degrees  of  freedom  and  many  coupled  excited  states 
such  calculations  may  be  prohibitively  complicated.   The  computational 
effort  required  for  a  coupled  channels  calculation  for  example 
increases  by  the  power  of  the  number  of  channels  involved. 


A  semiclassical  method  which  employs  classical  trajectories 
requires  much  less  computational  effort  than  an  exact  approach  for 
large  systems.  For  this  reason  it  is  important  to  develop  an 
efficient,  accurate  semiclassical  method  which  is  able  to  deal  with 
many  degrees  of  freedom  and  many  coupled  electronic  surfaces.   Of  the 
semiclassical  methods  mentioned  above  none  are  at  present  able  to 
effectively  deal  with  more  than  one  excited  state  potential  energy 
surface.  The  photodissociation  events  of  interest  here  have  a  bound 
initial  electronic  state  with  a  certain  set  of  vibrational  quantum 
numbers  and  a  final  electronic  state  with  a  set  of  final  vibrational 
quantum  numbers  for  the  molecular  fragment.   With  this  in  mind  another 
desirable  feature  of  any  photodissociation  theory  is  to  be  able  to  set 
initial  conditions  as  well  as  final  conditions  for  the  event,  i.e.  to 
calculate  cross  sections  for  state  to  state  transitions. 

The  formalism  in  the  following  treatment  is  designed,  then,  to 
include  certain  desirable  features.   These  features  can  be  summarized 
as  the  ability  to  efficiently  deal  with  (i)  many  degrees  of  freedom, 
(ii)  many  coupled  electronic  states,  and  (iii)  state-to-state 
transitions. 

The  particular  method  to  be  discussed  uses  the  self-consistent 
eikonal  approximation.   It  is  self-consistent  in  that  the  electronic 
transitions  are  self-consistently  coupled  to  nuclear  motions.  The 
eikonal  approximation  simply  means  that  the  nuclear  wave  function  is 
expressed  as  a  modulus  matrix  times  an  exponential  with  a  common 
eikonal  phase  and  that  the  short  wavelength  approximation  is  applied. 


The  general  process  of  interest  is  the  following  photodissociation 
reaction, 


M(v^)  +  *(k,(j)->  A(v^)  +  X(v^)  (I.l) 


where  M  is  a  polyatomic  molecule  in  a  given  electronic  and  nuclear 
state  \)j.,  A  is  an  atom  in  electronic  state  v.,  <)>  is  a  photon  with 
wavevector  k  and  polarization  a,   and  X  is  a  molecular  fragment  in 
electronic  and  nuclear  state  v^.   With  the  method  presented  in  the 
following  chapters  one  can  treat  photodissociation  events  as  in 
Eq.  (I.l)  involving  many  atoms  and  several  excited  potential  energy 
surfaces  in  a  way  which  is  computationally  practical. 


CHAPTER  II 
EVALUATION  OF  THE  PHOTODISSOCIATION  CROSS  SECTION 
USING  THE  EIKONAL  VAVEFUNCTION 


Radiation-Molecule  Interaction  in  the  Electron-Field  Representation 

The  interaction  of  a  photon  with  a  molecule  can  be  modeled  as  a 
collision  problem  [Micha  and  Svaminathan,  1985].   In  such  a  treatment 
we  consider  the  photon  to  be  occupying  a  finite  volume  V  in  space  and 
the  photon  field  to  be  contained  within  a  cube  of  side  length  L.   The 
molecule  M(Vu)  enters  this  space  and  interacts  with  the  field.   The 
field  is  quantized  so  that  the  molecule-field  system  is  described  by  a 
Harailtonian 

a=  H„  +  Hp  +  Hj^p  •  <^i-i> 

The  wavef unction  for  the  molecular  system  is  represented  by  |^M(Q)'*f 

where  the  "ket"  indicates  electronic  coordinates  and  Q  is  a  vector 

whose  components  are  mass  weighted  Cartesian  coordinates  for  nuclear 

motion.   The  Schrodinger  equation  for  the  molecular  Hamiltonian  H^  is 

n 

H„  |Y„(Q)>  =  E  \\(Q)>      ,  (II. 2) 

where  E  is  the  total  of  the  potential  and  kinetic  energy  of  the 
isolated  molecular  system.  In  all  subsequent  notation  the  use  of  the 
symbol  -  to  underline  a  character  will  represent  either  a  matrix  or  a 
collection  of  variables  whose  members  represent  components  of  a 
vector.  The  meaning  in  either  case  will  be  specified  or  obvious  from 
the  context. 
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The  second  term  in  the  Hamiltonian  H„  refers  to  the  radiation 

F 

field.   It  is  given  by  [Loudon,  1973] 


where  a,    and  a,    are  the  creation  and  annihilation  operators  of  a 

photon  with  wavevector  k,  polarization  a   and  energy  Hox   .   The 

product  of  the  annihilation  and  creation  operators  is  the  number 

operator  which  operates  on  an  eigenstate  |(kff)  )  of  the  radiation 

field  to  give  the  number  of  photons  in  state  (ka),  where  N  is  the 

number  of  photons.  The  photodissociation  process  treated  here  will 

involve  one  photon  present  in  the  initial  state  and  none  present  in 

the  final  state. 

The  Hamiltonian  term  H„_  refers  to  the  molecule-field 

nr 

interaction.   In  the  dipole  approximation  this  interaction  term  is 


^MF 


(r^)  =  -  J  D„(r'-r^)  .  E(r')  d^r'  (II. 4) 


where  r  is  the  position  of  the  center  of  mass  of  the  molecule  M.   The 
factor  Dj.  is  the  dipole  operator,  which,  in  the  center  of  mass  frame 
at  position  r  is 

5h(£>  =  S  ^a  Ea  *<£a  "  0  (II-5) 


where  {c  }  is  a  collection  of  charges  at  positions  {r  ).   The  electric 
field  operator  E(r)  in  second  quantized  form  is 


E(r)  =  E  g(")  e^[a    Xi^(r)  -  aj   Xj^Cr)*]  (II. 6) 

g(a))  =  i  [Ma)/(2So)]   ^^^,  (II. 7) 


where  So  is  the  vacuum  dielectric  constant  and  e  is  the  polarization 
unit  vector  with  a=x,y  for  linearly  polarized  light  and  ?=+,-  for 
circularly  polarized  light.  The  factor  X  in  Eq.  (II. 6)  above  is  the 
plane  wave  of  the  photon  inside  the  radiation  volume  and  has  the  form 


Xi^(r)  =  V  ^^^  exp(ik.r)   .  (II. 8) 


The  wavevector  k  has  components  2iijr/L  with  jr=0,l,2...  and  ^=x,y,z 
corresponding  to  cyclic  boundary  conditions  for  a  cube  centered  at  the 
origin  of  coordinates.   Replacing  Eqs.  (II. 5)  to  (II. 8)  in  Eq.  (II. 4) 
we  obtain 


-5A(^~a'!S)*V-^/2^"~'~^a^^^]        (II. 9) 
Di:,(£a.!s)  =  E  c^  £a  exp(ik.r^)  (11.10) 


The  Hamiltonian  in  Eq.  (II. 1)  acts  on  a  state  |j)  which  describes 
the  state  of  both  the  molecule  and  the  photon  field.  The  initial 
state  |i)  and  final  state  |f)  for  a  single  photon  event  are  given  by 


|i>  =  le^  M  (k  ff)>  and   |f>  =  Ie^  Ej.  a>  (11.11) 

where  p  and  p'  are  the  initial  and  final  momenta  of  the  center  of 
mass,  u  indicates  the  initial  molecular  state,  and  a  indicates  the 
final  state  of  the  fragments.  In  Eq.  (11.11)  the  state  |i> 
corresponds  to  a  molecule  initially  in  the  presence  of  a  photon  field 
and  the  state  |f>  to  a  dissociative  excited  state  with  no  photon 
present.  The  fragments  in  state  |f>  have  a  relative  final  momentum 
vector  of  gr-.   The  corresponding  total  energies  are 


E 


.  =  p  /(2m  )  +  E  +  Mw  (11.12) 


Ef  =  (E;)^/(2m^)  +  Efi/(2m)  +  E^  (11.13) 


where  m  is  the  total  molecular  mass  and  m  the  reduced  mass  for  the 
c 

relative  motion  of  A  and  X  of  Eq.  (I.l). 

The  dissociation  cross  section  is  the  quantity  of  interest  in  a 
photodissociation  experiment.  The  initial  and  final  states  are 
characterized  by  both  electronic  and  vibro-rotational  energies  of  the 
fragments.  Thus  the  cross  section  to  be  calculated  is  a  state-to- 
state  (electronic,  vibrational,  rotational)  resolved  quantity  at  a 
particular  photon  energy.  The  expression  for  this  cross  section  is 
[Rodberg  &  Thaler,  1967] 


da,. 


=  h-'   I  "^.c 


dS      i   J   «ic  . 


dPf.  Pf.  Rf.    ,  (II. lA) 


where  R,.  is  the  transition  rate  and  J.  is  the  incoming  photon  flux  in 
the  laboratory  frame  given  by 


J^  =  c/(2iiM)-^   ,  (11.15) 


where  c  is  the  speed  of  light.  To  obtain  Eq.  (11.15)  we  have  used 
momentum  normalization  for  which 


<r^|E^>  =  (2iiM)  2^^  exp  (ig^  .  r^/M)  .  (11.16) 


The  transition  rate  is  given  in  terms  of  the  transition  amplitude  T^. 
and  Dirac's  delta  function  with  respect  to  energy, 

^fi  =  r  '"^fil^  ^<^f  -  ^i^      •  (^^•^^> 

The  transition  amplitude  T^.  is  an  integral  over  electronic  and 
nuclear  coordinates  and  contains  an  operator  T  which,  for  low 
intensity  incident  radiation,  is  simply  the  molecule-field  interaction 
of  Eq.  (II. 9).   For  the  single  photon  process  to  be  considered  here 
only  the  first  term  in  H„p  is  appropriate.  A  single  photon  event 
corresponds  to  a  weak  electric  field  (i.e.  low  intensity).  The  form 
of  the  transition  amplitude,  then,  is 


Tfi  =  <f|T|i> 

=  -g(")  e^  •  <Efi°'|D,:,(r^,k)|w>  <E;|V  ^  e  ~'~^|e^>  .   (11.18) 


The  last  factor  containes  the  integral  over  r  which  gives 


ik.  •  r 
|<B'  |V-^/2  g  ~'~<^|d  >|2  ^  (2iiH)-^  S(E'  -  0     -   Mk)  .    (11.19) 


Therefore  the  integral  in  Eq.  (11.14)  becomes  independent  of  g  with 
the  restriction  that  p  '  =  p  +  Mk.  The  cross  section  now  contains  an 
integral  over  g,;.  only, 


d<Tj. 

j^-  =  J-^  g(w)^  (2II/M)  (2iiH)-^ 

X  Jdpj.  p2.  |<B^.alD;,(r^,k).e^|p>|2  S(E^  -  E.),   (11.20) 

where  p^.  is  the  magnitude  of  gr-.  If  we  note  that  Mk  <<  g  it 
follows  that  E  '~  E  •  With  this  very  reasonable  approximation  and 
using  E,  =  E.  we  get  from  Eqs.  (11.12)  and  (11.13)  for  p,. 


''fi 

2ni     u    a 


(11.21) 


Transforming  the  integral  in  Eq.  (11.20)  by  noting  that 

dE^  =  p,.  dpr./m  and  using  the  definition  of  J.,  the  differential 

cross  section  becomes 

'  2mnp^.  ^ 


^•^fi    ,  ,2 

d2-  =  e^"> 


Mc 


<Rfi«l5M(£a5°>-^al^^> 


(11.22) 


kiv 
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The  factor  in  brackets  will  be  denoted  as  D,,  indicating  that  the 
transition  is  from  a  molecule  initially  in  state  |i>=|y>  to  fragments 
in  a  final  state  |f>=|gj:.  a>  in  the  asymptotic  region.  The  final 
state  is  designated  by  the  states  of  the  atom  A  and  of  the  fragment  X 
asymptotically  as  as(r.,\)„)  and  v^s(r„,v)  where  T   refers  to  an 
electronic  state  and  v  to  a  collection  of  vibrational  quantum  numbers. 
In  this  treatment  we  will  assume  that  the  initial  bound  state 
wavefunction  |y>  for  the  molecule  in  the  center  of  mass  frame  is 
given.   It  remains,  therefore,  to  obtain  the  final  state  wavefunction 
a>. 

The  vector  Q  in  general  contains  components  for  the  relative 
position  R  of  A  with  respect  to  X,  the  internal  degrees  of  freedom  Q' 
of  X,  and  the  center  of  mass  position  r  .   In  the  following  notation, 
however,  we  will  be  working  strictly  in  the  center  of  mass  frame  and 
the  symbol  Q  will  refer  only  to  relative  and  internal  components. 

The  final  state  wavefunction  in  the  center  of  mass  frame  will  be 
expressed  in  the  coordinate  representation,  so  that 

<Qlefi  «>  =  I^M^^'Bfi)>     •  ("•") 

This  molecular  eigenstate  can  subsequently  be  expanded  in  a  basis  of 
the  electronic  states  *_  of  M  involved  in  the  process,  so  that 

\\(^)>   =  I!m(Q)>  iM^Q)  (11.24) 

|$„  >  =  (  |*->  |*„>  ...  |*,>  )  is  the  matrix  of  electronic 
wavefunctions,  and 


11 


14m 


\^ 


\ 


(11.25) 


is  the  matrix  of  nuclear  wavefunctions  where  h  is  the  number  of  the 

highest  state  energetically  accessible  in  the  reaction. 

We  will  express  Hj.  as  a  sum  of  the  nuclear  kinetic  energy 

operator  K   and  a  term  H^  which  contains  the  nuclear  repulsion  and 
nu  Q 

all  the  energies  of  the  electrons,  including  spin-orbit  coupling, 


M    nu    Q 


(11.26) 


K^^  =  (2M)-^  (-iHV)2 


(11.27) 


where  M  is  an  arbitrary  mass  which  comes  from  mass  weighting,  and  the 
symbol  7  is  a  multidimensional  gradient  whose  components  are  3/9Q., 
i=0  to  N.  The  gradient  7  to  the  second  power  is  defined  to  be 


2     2   2 
(7)^  =  E3^/3Q^ 


(11.28) 


The  coefficients  of  the  expansion  contained  in  vj;  are  the  nuclear 
motion  wavefunctions  on  the  electronic  potential  energy  surfaces. 
These  surfaces  are  defined  by  the  diagonal  elements  of 


Hq  =  <!mI  Hq  I!m> 


(11.29) 
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When  Eq.  (11.23)  is  substituted  into  Eq.  (II. 2)  and  left 
multiplied  by  <*„!  we  have  the  following  matrix  differential  equation 


2M 


'  ^  Iq  -  G(Q)  ^' 


^  §0  -  E  j  4„(Q)  =  0   . 


(11.30) 


We  will  Vork  in  the  diabatic  picture  [Smith,  1969]  in  which 


G(Q)=<!M(Q)l(>^/i)  ^I!m(Q)>  =  0 


(11.31) 


This  leaves  us  with  the  general  matrix  equation: 


1^ 
2M 


1 


^  %  -  eL(Q)  =  0 


(11.32) 


When  the  coupling  between  the  electronic  ground  state  (1=1)  and 
the  other  electronic  states  is  negligible,  H_  can  be  separated  by 
blocks, 


H  = 


^11 

0 

0   ....' 

0 

H22 

H23... 

\ 

0' 

0 

9 

• 

^32 

H33... 

D 

• 

• 
• 

•                        J 

(11.33) 


This  gives  us  one  equation  for  the  ground  state  wavefunction  and  a  set 
of  coupled  equations  for  the  excited  state  wavefunction.  Following 
this  argument  we  will  also  separate  the  wavefunction  expansion  as 
follows: 
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|T„(Q)>  =    IV  ti  .    E^  IV  *r 


(11.34) 


We  will  define 


\\r^^)>   =  l*i(Q)>  '^'i(Q) 


(11.35) 


to  be  the  ground  state  wavefunction,  and 


\^ex^9.)>   =  l!<Q)>  4(Q)  =  2  \^T^9.^>   +r^^^ 


(11.36) 


to  be  the  excited  state  wavefunction.   By  dividing  the  space  of  the 
electronic  basis  in  this  way  we  can  isolate  that  part  of  the 
wavefunction  which  involves  the  excited  states.   Thus  we  need  only 
solve  that  part  of  the  matrix  equation  which  deals  with  the  function 
1^  (Q)>'   The  matrix  equation  for  the  nuclear  wavefunctions  on  the 
excited  electronic  surfaces  is 


1^ 

2M 


^7 
1  ~ 


+  V  -  ieU  =  0 


(11.37) 


This  is  the  equation  we  solve  next  using  the  common  eikonal 
approximation. 


14 
The  Eikonal  Vavefunction 

The  nuclear  wavef unction  with  incoming  wave  boundary  conditions 
can  always  be  written  in  the  form  of  a  modulus  matrix  X(Q)  times  an 
exponential  with  a  common  phase  S(Q). 


!l;(Q)^  ^  =  ^(Q)  e^P[i  S(Q)/M]  (11.38) 


where  the  minus  sign  on  ijj  indicates  that  the  boundary  conditions 
chosen  are  those  of  an  incoming  wave.   We  now  substitute  Eq.  (11.38) 
into  Eq.  (11.37),  carry  out  the  T"   operation.  We  then  multiply  by 
exp[-iS(Q)/M]  to  obtain 


iM   2      iM  M^  2    (7S)^ 


2M 


(H  S))$  -  -^(VS  •  V)^)  -  2m2  ^  ^  -2M-^  "  <^  -  ^>^  =  0-  ^^^-39) 


We  can  also  write  this  as 


0  X  =  0  ,  (11.40) 


where  0  is  the  operator 


0  =  1/2M  (M/i  7  +  VS)^  +  V  -  E   .  (11.41) 


Note  that  0  is  in  the  form  of  a  Hamiltonian  minus  energy 


15 


0  =  H  -  E  (11.42) 

H  =  T  +  V  (11.43) 

T  =  T(p,7S)   ,  (11.44) 

where  the  normal  momentum  operator  p=M/i  7  has  been  replaced  by  the 

.sum  ^  +  7S  thus  defining  a  new  momentum  operator  p'  in  the  kinetic 
energy  operator  T, 


T  =(p')^/2M  .  (11.45) 


If  we  define  the  multidimensional  classical  momentum  to  be 


P  =  3S/3Q  5  7S  (11.46) 


we  get  the  following  for  Eq.  (11.39): 

0^  =  -2M  H^X  ^  f  ^M— ^  ^  r  15 ^  2M-  ^  ^  ^^  -  ^>^-    (^^-^^^ 

In  order  to  express  Eq.  (11.47)  in  the  form  of  a  Hamilton-Jacobi 
equation  we  can  carry  out  the  following  operations: 

>$""  0  X  +  (X^  0  X)^  =  0  (11.48) 


)C  0  )i  +  x"  0*  X=  0   .  (11.49) 


With  the  proper  manipulation  Eq.  (11.49)  becomes 
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(    1  ^ 

M 


2V 


X  +  fj^  P-tX^^  -  (?X)^X]  +  ^  X^X  -  2EX^)$  =  0.  (11.50) 


After  dividing  by  2)^  )^  this  can  be  put  in  the  form  of  the  Hamilton- 
Jacobi  equation, 


(P)  /2M  +7   =  E 
^  '•'  qu 


(11.51) 


where 


V  =  ^i*v-Y^-h  ^  ^  m  ^  m  ^ 


III-  mr  I 


.  (11.52) 


It  remains  to  solve  for  X(Q)'   To  do  this  we  first  substitute  Eq. 
(11.51)  into  Eq.  (11.47)  to  give 


%     l-ll-f-     P-?)$ 


^  ^X  +  (V  -  1V^^)X  =  0     (11.53) 


We  can  simplify  Eq.  (11.53)  by  introducing  the  following  notation: 


U  =  M/2M  V.P  , 


(11.54) 


Kj.,j.  =  5j.,j.(-H^/2M)7^Xj./Xj. 


(11.55) 


Equation  (11.53)  becomes 


•j,  P.7 

vV^+(V-l  V  )  +  K-  l(iU) 
1  M     ~   ~   qu    ~   ~^   ^ 


X  =  0 


(11.56) 
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Equation  (11.52)  becomes  the  expression  for  an  average  potential 
7  when  we  assume  the  second  term  to  be  negligible  and  make  the  short 
wavelength  approximation,  for  which  7'^)$=0  (see  appendix  A),  leaving  us 
with 


+  ...-1 


^   =  (X^>^)   ^*Y)$ 


(11.57) 


We  will  now  apply  the  short  wavelength  approximation  in  which  K  is 
negligible  and  V   becomes  V  to  give 


H 


P. 7 


■^  -—  +   (V  -   IV)  -  l(iU) 

in     ~   —     — 


)$=  0 


(11.58) 


In  chapter  III  we  will  solve  Eq.  (11.58)  using  trajectory 
equations  after  introducing  a  time  variable  within  the  new  set  of 
variables  A. 


CHAPTER  III 
INTRODUCTION  OF  A  TIME  VARIABLE:  TRAJECTORY  EQUATIONS 


A  New  Set  of  Variables 

The  matrix  equation  which  we  just  derived,  Eq.  (11.58),  will  be 
solved  numerically  using  classical  trajectories.   In  order  to  uniquely 
define  any  point  along  a  trajectory,  however,  we  need  to  introduce  a 
time  variable  and  indicate  the  values  of  each  component  of  Q  at  the 
initial  time.   Initial  values  of  the  Q  are  denoted  as  ^  and  are  given 
by  the  set 

g,  =  {q^^.  (\2'  '"  %^  (iii-i) 

where  q^  sets  the  initial  value  of  Q,  for  a  given  trajectory,  q„  the 
initial  value  of  Q„,  etc.  One  of  the  degrees  of  freedom  in  Q,  namely 
Q„,  will  be  set  by  the  relative  distance  between  A  and  X  which  will  be 
chosen  to  always  have  the  same  initial  value  for  each  trajectory.  The 
subscript  N  is  equal  to  one  less  than  the  number  of  degrees  of  freedom 
in  the  system  for  a  center  of  mass  reference  frame,  thus  there  will  be 
a  total  of  N+1  degrees  of  freedom  in  the  problem. 

Introducing  this  collection  of  variables  is  equivalent  to 
transforming  from  the  set  of  spatial  variables  Q  to  a  new  set  of 
variables  {t;^}  such  that  Q=Q(t;g).  The  set  g  together  with  the  time 
variable  will  be  denoted  A, 
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A  a  {t,  Qj^,  q2.  •••  qjq}  =  {t;  <i}   .  (III. 2) 

Each  trajectory  evolves  with  the  time  variable  t  separately.  The 
other  members  of  A  will  insure  that  each  trajectory  follows  a  unique 
path  in  Q-space  (even  though  a  particular  trajectory  will  cross  others 
many  times). 

It  is  convenient  to  define  sets  Q'  and  c^'  to  be  the  subsets  of  Q 
and  ^  respectively  which  refer  to  internal  variables.  They  will  each 
contain  N'  components 

Q'  =  {Qi',  Q2'---Q^,)  (III. 3) 

a'  ^  {Qi'.  <l2'-"%>)  (III. A) 

where  N'  is  the  number  of  internal  degrees  of  freedom  in  the  system. 
We  can  write  the  initial  value  of  any  the  Q.'  in  terms  of  an  angle  q.' 
defined  by 


Q.'(t=0)  =  (y.)^^^  A^.  sin(q.')  (III. 5) 


where  A  .  is  the  maximum  value  of  Q.'  for  quantum  number  v.  and  y.  is 

VI  1^1        1 

the  reduced  mass  for  the  i   internal  mode.   Note  that  for  the  Q' 
corresponding  to  normal  vibrational  modes,  if  the  potential  energy 
along  this  Q.'  is  harmonic  in  the  asymptotic  region  (Qr-)^")  then  we  can 
write 


Q.'(t;^)  =  (y.)  ^^2  A^.  sin[(o.  ( t-t°)  +  q.']   ,       (III. 6) 
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where  w.  is  the  frequency  along  Q'.    and  t  is  the  initial  time.   The 
asymptotic  region  is  that  part  of  Q-space  where  the  potential  is 
separable  in  Q' ,  i.e. 

3V(Q)/3Q'  =  f(Q')  (III. 7) 

Thus  Eq.  (III. 6)  gives  us  the  functional  relationship  for  a 
vibrational  mode  between  Q'  and  c['  for  a  harmonic  potential  in  the 
asymptotic  region. 

Now  that  we  have  introduced  the  time  variable  we  would  like  a 
prescription  for  determining  the  trajectories.   The  classical 
equations  of  motion  for  position  and  momenta  can  be  formulated  using 
Hamilton's  equations.   We  define  Eq.  (11.51)  to  be  the  Hamiltonian  H 
and  impose  the  restriction  of  energy  conservation  to  obtain 


gf=0   .  (III. 8) 


This  leads  to  Hamilton's  classical  equations  of  motion 


^   31 
dt  "  3P 

dP 


(III. 9) 


^"^        dr  =  -lr  •  (i^i-io) 

Both  Q  and  P  can  be  written  as  a  collection  of  variables, 

Q  =  {Qq,  Qy   Q2---Qn)  (III. 11) 

P  =   {Pq,  P^,  P2...Pn^   •  (III. 12) 
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For  each  component  of  Q  and  P,  then, 


It  is  helpful  at  this  point  to  discuss  mass  weighted  coordinates, 
We  define  the  following  to  be  mass  weighted  coordinates, 

Q.  =  nr    X.  /  fir  dii.is) 

Pyi  ; 

p.=  — —  fir  (III. 16) 

where  X.  and  P„.  are  the  original  Cartesian  coordinates  for  the  i 
degree  of  freedom.   The  mass  M  can  be  arbitrarily  set  to  some  value 
without  affecting  the  results  of  the  derivation.   From  Eqs.  (III. 15) 
and  (III. 16)  it  follows  that 


dQ.(A) 
Pi(A)  =  M   -ji .  (III. 17) 


Equation  (III. 17)  shows  the  advantage  of  using  mass  weighted 
coordinates  in  that  the  reduced  mass  y.  does  not  appear  explicitly  in 
the  relationship  between  P  and  Q. 

Turning  now  to  the  differential  equation  derived  in  chapter  II  we 
can  simplify  Eq.  (11.58)  using  the  definition  of  V  in  chapter  II  and 
noting  that 
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P.V     N  dQ.   g      ^ 

IT  =  .^Q  dT  3q7  =  dt  •  (III. 18) 


The  Eikonal  Vavef unction  Along  a  Trajectory 
Having  introduced  the  set  A  we  can  rewrite  )^  using  this  new  set 


as 


-  I     \   Y(t')  dt' 


)$[Q(A)]  =  X  C(A)  e      t"  (III.  19) 


where 


Y  =  J  7.P  •  (III. 20) 


and  X  is  the  value  of  X  at  t  which  is  set  by  initial  conditions.   If 
we  replace  Eqs.  (III. 18)  and  (III. 19)  into  Eq.  (11.58)  we  get 


(M/i)  dC/dt  +  (V  -  1V)C  =  0   .  (III. 21) 

All  that  remains  is  to  solve  Eq.  (III. 21)  for  C  noting  that  the 
time  dependence  is  through  trajectories.  As  we  travel  along  a 
particular  trajectory  we  will  be  solving  for  C  and  V  simultaneously 
because  the  form  of  V  is  now 

V  =  (C'^C)~-'-C"^VC   .  (III. 22) 


23 

We  will  also  construct  the  nuclear  wavef unction  along  the  trajectory. 
It  now  has  the  form 


r~  (Q)  =  X(Q)  e       =  X  C(t;a)  e  ^  t°     e"       (III. 23) 


where  S(Q)-»<T(t  ;^) .  It  will  later  be  shown  that  the  quantity  involving 
Y  simplifies  to  a  function  of  the  momenta  and  other  terms  which  can  be 
generated  numerically  along  the  trajectory. 

The  final  part  of  ■jr~'   left  to  construct  is  the  factor 


e 

To  do  this  we  refer  to  the  definition  of  P  in  Eq.  (11.46)  and 
integrate  to  get 

TdS  =    rP  .  dQ  (III. 24) 


S(Q)  -  S(Q°)  =  E     I  Pi  ^Qi  (III. 25) 

i=OQoJ 


,  N     -,t  „ 
S[Q(t;^)]  =  S[Q(t°;^)]  +  i  ^      ?•  dt  h  a(t;^)      (III. 26) 

i=0  o 

where  Q  indicates  all  initial  values  of  Q  and  t  the  initial  value  of 
t.   Thus,  according  to  Eq.  (III. 26),  to  calculate  (i(t;^)  we  need  only 
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integrate  the  square  of  the  momenta  along  each  trajectory  over  the 
time  variable. 

We  will  now  develop  a  prescription  for  generating  y.      The 
function  y  =  H*~'''^  would  seem  to  pose  quite  a  problem  due  to  the  fact 
that  it  is  in  the  exponent  and  contains  partial  differentials.   It 
turns  out,  however,  that  the  entire  exponent  can  be  reduced  to  a  very 
simple  expression  containing  the  Jacobian  for  the  variable 
transformation  from  Q  to  {t;^}. 

Let  us  first  recall  that  the  term  involving  y  appeared  in  y^~'   as 
a  result  of  having  transformed  the  differential  equation  from  spatial 
coordinates  Q  to  the  set  of  coordinates  A  which  contains  a  time 
variable.  Using  the  definitions  of  P  and  7  the  product  in  y   becomes 


P(Q) 


N 

E 
i=0 


3P. 
: 

3Q, 


(III. 27) 


Let  us  look  at  one  term  in  the  sum  on  the  right  side  of  Eq.  (III. 27), 
i.e.  3P./3Q..  We  know  from  calculus  that  we  can  always  relate  a 
partial  derivative  in  one  coordinate  system  with  the  partials  in 
another  coordinate  system  providing  we  can  construct  the  Jacobian  of 
the  transformation  between  the  two  systems.   In  this  case  the 
particular  relation  between  partial  derivatives  is  the  following: 


3P. 

1 

3Q7 


J 

t      q^...q....q^ 

J 

}    ^r--'ii---V 

(III. 28) 
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Note  that  the  i   column  of  the  numerator  in  Eq.  (III. 28)  contains  P. 
rather  than  Q..  The  explicit  form  of  the  Jacobian  in  Eq  (III. 28)   is 


_t  ...  Qi-.-QN^ 


3Q0 
3t 


3t- 

!^ 

3qi 


3Q, 


9q 


N 


30 
3t 
3Q 
3^ 


N 


N 


3Q 
3q^ 


(III. 29) 


If  we  use  Eq.  (III. 28)  in  Eq.  (III. 27)  the  product  appears  as  the 
following  sum  over  i  of  the  determinant  in  Eq.  (III. 28)  divided  by  the 
Jacobian, 


N 


i=0 


3Q(- 

FT 

30, 

3qi 


0 


3Qf 
W 


N 


3P. 

1 

9t 

3P. 

i_ 

3qi 


3P, 


3q 


N 


30 
Ft 
30 
9qi 


N 


N 


3q 


N 


(III. 30) 


where  J  is  the  Jacobian  of  Eq.  (III. 29).   The  summation  of 
determinants  in  Eq.  (III. 30)  calls  to  mind  the  expression  one  obtains 
when  taking  the  derivative  of  a  single  determinant.  The  derivative  of 


the  Jacobian  determinant  with  respect  to  q .  is 
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N 

M_  _  V 
aq.  ^  L 

J    i=0 


3Q, 


3^0, 


3Qk 


3t      •• 

••    3q.3t      *••• 

3t 

3Qo 

.\ 

3Q 

8qi    •• 

"    Sq^aqi   '*•• 

3q 

3^0 

• 

2 
3^0    . 

3Q 

^%    " 

"    3^j3^N   ■'■' 

3q 

(III. 31) 


The  determinant  in  Eq.  (III. 31)  looks  very  similar  to  that  in  Eq. 
(III. 30)  except  that  the  i   column  contains  double  derivatives 
rather  than  the  derivatives  of  P..   If,  however,  we  recall  that 
P.=M(dQ./dt) ,  then  we  can  make  these  determinants  identical  by 
choosing  q.=t  to  obtain 


3J_ 
3t 


N 

z 

i=0 


3Qq 

^\ 

3Q 

3t      ** 

"        3t2      •• 

••    3t 

3Qo 

3^0, 

3Q 

3qi    •• 

••    3t    3q^    •• 

"    3q 

3Q0 

3^Q    . 

3Q 

3q 


N 


3t  3q, 


^q^ 


(III. 32) 


If  we  now  switch  the  order  of  differentiation  on  the  double 


derivatives  and  replace  3Q./3t  by  P./M  we  get  that 


N 

3t  -  m/_, 
1=0 


3Q0 
3q^ 


3Q(. 
3a 


N 
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9P. 

1 

3t 

3P. 

1 

3qT 


3P. 


3q 


N 


N 
3q^ 


(III. 33) 


3Q, 

3^ 


After  comparing  Eq.  (III. 33)  with  Eq.  (III. 30)  we  can  express  y  in 
terms  of  the  time  derivative  of  the  Jacobian  which  gives 


U^d^.hJ 


dt  '  ^H 


(III. 34) 


1  ^   ^_  1  dJ  _  d  ln(J) 


J  dt 


dt 


(III. 35) 


Recall  that  the  term  involving  y   appears  in  the  nuclear 
wavefunction  ;|;    as  an  integral  in  the  exponent.   Replacing  Eq. 
(III. 35)  in  the  integral  gives 


- 1  f:^ 


t')  dt' 


-2  f ^  1" 


(J) 


-\[    ln(J)  -  ln(J°)  ] 


expj   -  ^   [    ln[  gi  arg(J)]    +  In    |j  |   -  In  [^i   arg(J°))    -  ln|j°|] 


1    ,. 


expj   -J    [i   arg(J)   -   i   arg(J   )] 


1/2 


(III. 36) 
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where  J°  refers  to  the  value  of  the  Jacobian  at  time  t  .  Thus  we  have 


4  Jr(t') 


dt' 


1/2 


exp{-^[  arg(J)  -  arg(J°)  ]}.  (III. 37) 


Our  nuclear  wavefunction  can  now  be  constructed  along  a 
trajectory  defined  by  A  a  {t;^'}.  The  final  form  of  this  wavefunction 
is 


i^  \t;a)  =  X°  C(t;a) 


J"(0;a) 


J(t;a) 


1/2   1 


-2U(t;a)  R  a(t;a)   (111.38) 


,0  . 


where  J  is  the  Jacobian  at  time  t,  J  is  the  Jacobian  at  the  initial 
time,  and  y  contains  the  argument  of  the  Jacobian  as  follows: 


M(t)  =  arg(J)  -  arg(j") 


(III. 39) 


The  amplitude  matrix  C  is  self-consistently  solved  for  through 
Eq.  (III. 21).   The  factor  X  is  a  constant  whose  value  will  be 
determined  in  chapter  IV  by  initial  conditions.   The  other  terms  are 
constructed  from  quantities  generated  along  the  trajectories.  The 
nuclear  wavefunction  calculated  in  this  manner  will  be  used  in  the 
next  chapter  to  calculate  the  transition  integral  which  appears  in  the 
cross  section  of  Eq.  (11.22). 


CHAPTER  IV 
THE  TRANSITION  INTEGRAL  IN  A- SPACE 


Transformation  to  an  Integral  Over  the  Time  Variable 


The  quantity  of  interest  to  the  experimentalist  is  the  state-to- 
state  cross  section.   It  is  calculated  using  a  transition  integral 
between  the  ground  state  wave  function  v|/,  and  the  function  jj/^   .   This 
is  calculated  from  the  quantity  D,.  which  appears  in  Eq.  (11.22). 

The  general  form  of  the  transition  integral  is 


°fi  =  I  fif  ^^^5-^  ^i  ^'^  ^^^'^^ 

where  T  refers  to  integration  over  all  space.  The  subscript  f  has 
been  introduced  on  <J;    to  indicate  the  asymptotic  form  of  the 
wavefunction.  The  symbol  D  is  a  function  of  nuclear  variables  only. 
It  is  derived  from  the  integral  of  initial  and  final  electronic  wave 
functions  and  the  dipole  operator  D'  over  electronic  coordinates.   It 
has  components 


[5(Q)]r'i  =  <*r'^5^l5Ml*i^^^>-  ^^^'^^ 


The  integration  over  space  dx  in  Eq.  (IV. 1)  can  be  written  as  a 
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30 

multiple  integral  over  the  set  of  coordinates  Q  a  {Qq,  Q,,  ...  Q^) 
Thus  the  transition  amplitude  is  a  multiple  integral  in  (N  +  1) 
variables  which  can  be  written  as 


''fi  =  Jl4f"^l'5-S'^i 


dQ 


(IV.3) 


In  order  to  calculate  this  integral  using  trajectories  we  must 
again  make  a  variable  transformation  from  Q  to  A  as  in  the  last 
section.  The  integral  in  Eq.  (IV.3)  becomes 


°fi  =  J  lif"^]'  °---  h 


'%^i  •••V 


t    qi 


^Nj 


dA 


(IV. A) 


where 


^^0^1 


is  the  absolute  value  of  the  Jacobian  of 


transformation  and  the  collection  of  variables  A  is  defined  by  A 
=  {t,^}  where  ^={q-,  »qo> . .  .qjj}  contains  the  initial  coordinates  for  a 
trajectory.  To  show  explicitly  that  D^.  now  contains  an  integral  over 


time  we  write 


°fi  =  J  tif  ^J^5-e  *i  |J|  dt  d^ 


(IV. 5) 


We  can  approximate  the  integral  over  d^  by  a  sum  over  ^.   This  leaves 
us  with  a  sum  of  integrals  over  dt: 


°fi  =  V^^n)   j'lin'^l\-~  *i.n  l^n  I  ^^  ^  ^ /^fi.n  ^^^'^^ 
n=l        0  n=l 
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D^:.   =  (Aq  )    r  [ij;^"^]*D  .e  v|/.    |J  I  dt   .    (IV. 7) 


For  simplicity  we  have  dropped  the  subscript  f  from  \      •     The 

label  (n)  which  appears  as  a  subscript  in  the  summation  refers  to 

individual  trajectories.  Therefore,  ijj  '  is  a  function  which  is 

generated  along  the  n   trajectory.   Since  each  trajectory  is  defined 

by  its  values  of  q  it  can  be  said  that  (n)  labels  each  unique  set  q  . 
**  '*n 

This  can  be  seen  by  defining 


ii  ^O  a  i^'^tjq^^)   .  (IV. 8) 


According  to  Eq.  (IV. 6)  each  trajectory  n  will  contribute  a 

certain  amount  D,.    to  the  cross  section.  This  integral  will  be 
ri,n  ° 

calculated  along  the  n   trajectory  by  simultaneously  generating  each 
term  in  the  integrand  and  integrating  step  by  step  over  t.   In  order 
to  do  this  the  integrand  must  be  composed  of  terms  which  can 
themselves  be  generated  step  by  step  over  t. 


Generating  the  Jacobian  Along  a  Trajectory 


In  chapter  III  we  developed  the  detailed  form  of  ;|;   •   It  was 
also  shown  how  the  factors  C(t)  and  a(t)  could  be  generated  along  the 
trajectory.  The  other  factors  appearing  in  ^~'   contain  the  Jacobian 
as  an  argument.  This  is  the  same  Jacobian  as  that  which  appears  in 
the  transition  integral  of  Eq.  (IV. 4)  after  transforming  from  dQ  to 
dA.   It  has  not  yet  been  shown  how  J  can  be  calculated  along  each 
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trajectory.  This  will  be  done  next  as  we  discuss  a  way  of 
simultaneously  generating  all  of  the  terms  in  J. 

The  Jacobian  of  interest  is  the  same  as  that  in  Eq.  (III. 29). 
The  object  is  to  numerically  solve  for  each  term  and  then  to  carry  out 
the  operations  of  the  determinant.  This  is  to  be  done  at  each  time 
step  along  a  given  trajectory  in  order  to  obtain  the  Jacobian  for  that 
trajectory  as  a  function  of  time. 

To  generate  numerically  the  terms  in  Eq.  (III. 29)  we  need  to 
somehow  express  the  time  derivative  of  each  term  as  a  function  of  the 
other  terms  or  of  known  functions.   First  let  us  define  the  quantities 


Q^  a  9Q./9q.    and   P^  =   3P./8q. 


(IV. 9) 


where  j  =  0,  1,  2,    ...N;   i=  1,  2,    . . .N.   Using  the  first  definition 
in  Eq.  (IV. 9)  the  Jacobian  now  appears  as 


%  Qr 


t  q 


•V 


1-"4nJ 


1 
M 


N 

1 
N 


(IV. 10) 


To  obtain  the  time  derivatives  of  P.  and  Q.  let  us  start  by  recalling 
Eqs.  (III. 13)  and  (III. 14).   If  we  take  the  derivative  of  the  left 
side  of  Eq.  (III. 13)  with  respect  to  q.  and  switch  the  order  of 
differentiation  we  get 
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2       2       i 
3  Q.     3  Q.     3Q. 

±_     =  l_   =   ^       .  (IV. 11) 

3q.3t    3t9q.     3t 

Taking  the  derivative  of  the  right  side  of  Eq.  (IV. 11)  with  respect  to 
q.  gives  the  relation 

3Q^/3t  =  3P./3q.  =  P^   .  (IV. 12) 

Thus  Eq.  (IV. 12)  gives  us  one  of  the  time  derivatives  we  need.   If  we 
similarly  take  the  derivative  of  the  left  side  of  Eq.  (III. 14)  with 
respect  to  q.  we  get 


2       2        i 
3  P.     3  P.     3P. 

l_     =  ^_     =   ^     .  (IV. 13) 


3q.3t    3t3q, 


3t 


Following  the  same  procedure  as  before,  we  set  Eq.  (IV. 13)  equal  to 
the  derivative  of  the  right  side  of  Eq.  (III.IA)  as  follows: 


'^]  3^  V 

-IT    =  -  3^7307    •  <^^-^^> 

However,  since  we  have  no  expression  for  the  derivative  of  7  with 
respect  to  ^,  an  additional  step  must  be  taken.   Ue  know  the  quantity 
7  only  as  an  explicit  function  of  Q.   We  also  know  that  the  derivative 
dependence  of  Q.  on  q.  is  one  of  the  terms  being  generated 
(Q^=3Q. /3q.).  Therefore  we  can  expand  3  V/3q.3Q.  in  terms  of  the  Q^ 
by  using  the  chain  rule. 
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It  is  important  to  note  that  the  average  potential  7  is  a 
function  of  the  components  C_(t;^)  of  the  amplitude  matrix  C(t;^)  from 
Eq.  (III. 22).   It  is  convenient  to  express  C-,  in  terms  of  its  real  and 
imaginary  parts 


Cj.  a  K  (Xj,  +  iYj.) 


(IV. 15) 


where  K  is  a  constant  to  be  specified  later.  With  the  definition  of 
Eq.  (IV. 15)  we  note  that 


7  a  7(Q;  X,Y) 


(IV. 16) 


where  X  and  Y  are  sets  containing,  respectively,  the  collection  of  X_ 
and  Y^  components  of  all  C_.   Since  V  is  a  function  of  the  Xy,  and  Y_ 
these  must  be  included  in  the  chain  rule  expansion. 

Define,   F.  s   -37/9Q.  (the  force  along  Q.)  and  substitute  into 
Eq.  (IV. 14).   Using  the  chain  rule  we  get 


3F. 


3q, 


N 

E 

k=0 


9F.^ 


y^Qj 


3Q, 


k'C 


9q,      r 


J 


^dX— y 


3X, 


Q,Y3^i 


raF.^ 


l9Yj,; 


9Y, 


Q,Y3^i 


(IV. 17) 


In  Eq.  (IV. 17)  we  have  introduced  two  new  partial  differentials  which 
will  also  be  numerically  generated  along  each  trajectory.   We  will 
give  them  the  following  notation: 


9X, 

X^  =  ^ 
^r  -  9q, 


9Y_ 

y'  =  -£ 
r   9q. 


(IV. 18) 
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Using  Eqs.  (IV. 17)  and  (IV. 18)  and  the  definition  of  F.  we  get  the 
following  for  the  time  derivative  of  P.: 


3P-: 


at 


N 

E 

k=0 


3^V 


3Q^3Qj^ 


l9_ 
3X, 


3V 


3Q 


V 


x^  +  ^ 
r  *  3Y, 

X,Y 


''3V 


3Q, 


(IV. 19) 


We  have  now  expanded  the  partial  derivative  of  Eq.  (IV. 14)  in  terms  of 

quantities  which  we  know  how  to  calculate;  however,  we  had  to 

introduce  the  factors  X-,  and  Y-,  .  Of  course,  we  have  no  analytical 

expression  for  these  factors  so  we  must  again  develop  expressions  for 

the  time  derivative  of  each  and  generate  them  numerically  as  we  did 

with  Q.  and  P"^. 
J      J 

The  first  step  is  again  to  look  at  the  time  derivative  of  the 
function  whose  partial  derivative  we  are  considering.   In  this  case  we 
need  the  time  derivatives  of  X_  and  Y_.   To  get  an  expression  for 
these  time  derivatives  we  use  Eq.  (III. 21)  and  separate  the  real  and 
imaginary  parts  given  in  Eq.  (IV. 15)  to  obtain 


dX, 


3V 


dt  -  ^T  -  ^   ^TT'    ^V    ~   3Yj. 


(IV. 20) 


dY 


I  =  Y  -  -V  U    X   -  ^^ 

dt   "  r  "  ^      TV      T'    ~    3Y^ 


(IV. 21) 


where  the  matrix  W  is  defined  as  W=V-1V  and  we  assume  the  matrix  V 
defined  in  Eq.  (11.33)  to  be  real.  Next  we  take  the  partial 
derivative  of  Eqs.  (IV. 20)  and  (IV. 21)  with  respect  to  q.  and  switch 


the  order  of  differentiation 
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3X, 


3q.    3q. 


Y 

rr'  r' 


dxj 

dT 


s  x^ 


(IV. 22) 


raw 


rr' 


[3qi 


"•■  rr'  F' 


(IV. 23) 


Once  again  we  run  into  a  situation  where  we  must  expand  a  partial 
derivative  with  respect  to  q..   By  using  the  chain  rule  expansion  of 


9W 


rr' 


3q. 


we  obtain 


Xt  = 


T'      M  k 


■*  rr' 


3Q,, 


oi* 


x,y 


r" 


raw       9w 
"'^rr'  i   ""rr'  i 

3x^..  ^r"  *   3Y^..   r" 


T" 


T" 


+  2j  ^yv  ^T' 


(IV. 24) 


The  expression  in  the  equation  above  can  be  simplified  by  noting  that 
the  only  terms  in  W  that  are  dependent  on  X  and  Y  are  those  containing 
V,  namely,  the  diagonal  terms.   We  can  indicate  this  by  using  a  delta 
function  S__,  to  get  the  final  expression 


x;  = 


E    E 
r'  I  k 


raw 


rr' 


9Q, 


Q,.  - 


hi 


^TV   E 


3V 
3X 


r" 


T" 


3V 
3Y 


r" 


T" 


"^  E  ^rr'  ^r' 


(IV. 25) 
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At  last  we  have  an  expression  for  which  all  terms  can  be  calculated 
(either  analytically  or  numerically).   Using  similar  steps  for  Y-  we 
arrive  at  the  time  derivative 


-Yt  = 


E  E 

r'  I  k 


3W 


rr' 


3Q,, 


^■nr^r    Ll 


TF' 


^X,Y 


r" 


9V_  yi   ^  3V_  i 

axp,  ^r"  *  9Yp/r" 


K> 


+  2]  Wpj.,  Xj., 


(IV. 26) 


Returning  to  Eq.  (IV. 19)  we  see  that  having  derived  Eqs.  (IV. 25) 

and  (IV. 26)  we  have  an  expression  for  the  time  derivative  of  P.  and 

thus  for  Q.  in  terms  of  quantities  we  can  calculate.   We  are  already 

generating  all  P.  in  time  and  so  we  have  all  the  terms  to  calculate  J 

along  the  trajectory.  The  only  terms  in  D^.,   of  Eq.  (IV. 7)  which  we 

have  not  yet  talked  about  are  those  which  we  need  to  be  given,  namely, 

D  .e  and  ti 
-n  ~     l,n 

The  ground  state  wave  function  <j;.  will  be  given  as  a  function  of 
Q;  therefore,  ^.  (t)  is  simply  the  value  of  ^.  at  that  point  Q  which 
corresponds  to  the  n   trajectory  at  time  t;  i.e. 


^.n  =  ntQ~(^'VJ 


(IV. 27) 


The  factor  D  .e  will  similarly  be  given  as  a  function  of  Q. 

~n   ~  J  a 

We  now  have  all  the  elements  required  for  the  transition 
integral.  The  entire  expression  in  detailed  form  is 


D,.  =  E  D.. 
fi   n  11, n 


(IV. 28) 
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where 


D,.   =  (/iq  )   X^  C*(t) 
fi,n  ^^^n'  J-'ti  ~n^ 


n 


J 
n 


^''^     i[arg(J^)  -  arg(J°)]/2  -ia^(t)/M 
e  e 


X  g((o)D^(t).e  <^^^^i^)      \\\     dt  (IV. 29) 


The  symbol  J  is  defined  as  the  value  of  the  Jacobian  at  time  t   (the 
initial  time)  for  the  n   trajectory. 

Notice  in  Eq.  (IV. 29)  that  there  will  be  a  singularity  in  '^~' 

1  /  0 

caused  by  the  term  |j  /J |    each  time  J  passes  through  zero. 
Remember  that  ^~     is  constructed  using  the  short  wavelength 
approximation.   This  approximation  breaks  down  whenever  P  goes  to 
zero.   This  happens  as  a  trajectory  approaches  a  caustic  where  J  goes 
to  zero  and  »j;    passes  through  a  singularity.   Fortunately  these 
singularities  do  not  appear  inside  the  transition  integral  of  Eq. 
(IV. 29)  because  of  the  Jacobian  which  comes  from  tranforming  from  dQ 
to  dA.  The  fact  that  the  singularities  in  ^~     are  integrated  out  in 
this  way  makes  it  numerically  possible  to  use  the  semiclassical  wave 
function  t    without  having  to  integrate  through  poles. 


Asymptotic  Conditions 


There  is  still  a  problem  in  the  way  ij;    behaves  asymptotically. 
In  order  to  satisfy  asymptotic  conditions  the  nuclear  wave  function 
must  have  the  following  form  as  R-*«°  and  t-»t  : 
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1  p°.R 

^^PiQ)   =  (2ixM)"^^^  u^(Q')  e^  ~^'~  (IV. 30) 

where  d  is  the  the  dimensionality  of  the  space  that  A-X  is  constrained 
to  (d=l,  2,   or  3),  u  (Q' )  is  the  internal  nuclear  motion  wave  function 
(including  internal  rotation)  of  the  fragment  X  in  a 
vibrational/rotational  state  with  quantum  numbers  v,  Q'  is  the  subset 
of  Q  which  contains  all  the  internal  modes  of  the  fragment  X,  and  R  is 
the  vector  corresponding  to  the  relative  displacement  of  atom  A  from 
the  center  of  mass  of  the  fragment  X.   To  get  f"~     asymptotically  we 
will  look  at  how  each  term  in  Eq.(IV.29)  behaves  as  Q-Kl  ,  where  Q 
refers  to  values  of  Q  which  are  in  the  asymptotic  region. 

Asymptotically  we  can  express  the  eikonal  wave  function  as 


1^-'  -  x°  % 


a 

J 
a 


1/2  4  y„  la 


^     ^  e^  ^  (IV. 31) 


where  the  subscript  "a"  indicates  that  the  various  functions  are  in 
the  asymptotic  region.   To  solve  for  X  C  let  us  recall  Eq.  (III. 21) 
and  look  at  the  term  (V  -  IV).  Asymptotically,  the  off  diagonal 
elements  of  V  will  be  zero.   If  we  choose  the  wave  function  initially 
to  be  entirely  on  excited  surface  r,  then  the  term  V   becomes 


V  =  V  =  V™  s  E^  .  (IV. 32) 

Q-)Q  ^    ^^    ^ 
~  ~a 
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The  elements  in  V  -  17  =  W  then  become 


[W  ]. .  =  5. .(E.  -  E_)  (IV. 33) 


thus  having  off  diagonal  coupling  elements  equal  to  zero  and  a  zero 
for  the  r   diagonal  element. 

Because  there  are  no  off  diagonal  elements  in  the  matrix  equation 
asymptotically  and  because  C  will  have  the  following  form 


C 
-a 


?r 
lo  J 


a  C  -  ,  (IV. 3A) 

d.,  i 


Eq.  (III. 21)  will  reduce  to 


dCj. 
?    dT  =  °  '  <^^-35) 

thus  C  T"  is  constant  in  time. 

The  Jacobian  has  a  very  interesting  and  simple  form 
asymptotically.   In  the  asymptotic  region  any  internal  motion  along  Q. 
will  not  affect  the  relative  motion  along  Q.  and  vice  versa.   Using 
this  fact  we  can  set  the  terms  3Q./9q'  and  3Q'. /3q.  equal  to  zero  where 
Q.  and  q.  refer  to  relative  variables.  The  Jacobian  becomes 
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1 


0 
3Q 


0 


8Q, 


8Q 


1 


9qi   3qi 


8Qi 


3q„    3q, 


2 
3Q. 

9q] 
3Q. 
3q: 


0  . 

0  . 
3Q- 

3q^ 


0 

0 

3Q, 
3? 


N 


(IV. 36) 


3Q- 

3a" 


N 


'% 

3qN 


with  the  Jacobian  determinant  having  a  block  for  relative  variables 
(i=0,l,2)  and  a  block,  for  internal  variables.  Thus  the  determinant 
reduces  to  a  product  of  the  Jacobians 


J  =  J  „ 
a    a,R 


(%   Ql  ^2' 
^1  ^2J 


^1  ^2  ■■■  ^N' 


l^i^2 


%') 


(IV. 37) 


Outside  the  interaction  region  the  potential  along  Q^,    Q. ,  and  Q»  does 

not  change  and  so  P^,  P,,  and  P^  are  constant.   In  all  of  the 

following  equations  a  superscript  o  indicates  a  value  at  the  initial 

time.   If  we  restrict  the  definition  of  q,  and  q„  such  that  q,  =  Q, 

and  q„  =  Q^,  then  the  determinant  symbolized  by  J  „  becomes  simply 
L  L  a  y  K 

the  constant  P^^.   An  analysis  of  the  Jacobian  J'  for  the  asymptotic 
internal  degrees  of  freedom  shows  that  it  will  be  dependent  on  the 
combined  frequencies  of  the  normal  modes  of  vibration  for  the 
fragment.  A  specific  example  of  the  asymptotic  Jacobian  will  be 
presented  in  chapter  V  for  the  simple  case  of  a  system  with  two 
degrees  of  freedom.   For  the  present  discussion  it  is  sufficient  to 


note  that  J  is  a  periodic  function  dependent  on  an  effective 
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frequency  w  which  is  itself  dependent  on  the  frequencies  of  the  normal 

modes.   If  T  is  the  period  of  the  internal  modes  then  the  value  of  J 

a 

at  t  +T  is  equal  to  its  value  at  t  . 


Based  on  the  discussion  above,  the  asymptotic  Jacobian  has  the 


form 


V^'^>  =  ^0  ^4<^'^> 


(IV. 38) 


where  the  (t;c[)  dependence  of  the  Jacobian  has  been  indicated.   With 
this  Jacobian  the  asymptotic  nuclear  wave  function  has  the  form 


,.(-) 


=  * 


(-) 

r,a 


x°  c, 


J:(t";a') 


J:(t;a') 


2       -J   U(t;^)   jj  <T(t;q) 


(IV. 39) 


The  next  term  to  look  at  is  the  action  S.  To  see  how  S  behaves 
asymptotically  we  can  rewrite  Eq.  (III. 24)  as 


S(Q)  =  S(Q°)  +   ['  P  .  dQ 


(IV. 40) 


The  momentum  P  is  a  multidimensional  vector  which  can  be  separated 
into  the  sum  of  two  vectors, 


P  =  Pj,  +  P' 


(IV. 41) 


where,  P„  refers  to  the  relative  momentum  between  the  atom  A  and 
fragment  X,  and  P'  to  the  internal  momentum  of  the  fragment  X.   With 
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these  definitions,  then,  the  action  integral  becomes 


Q)  =  S(Q°)  +  rPj^  .  dR  +   ["P'  .  dQ'  (IV. 42) 


S(Q) 


This  can  be  written  as  follows  to  isolate  the  term  P«  .  R, 


S(Q)  =  P°.R  -  Pr.R°+  S(Q°)  +  "(Pj^-  P°).dR  +  ~  P'.dQ'  .(IV. 43) 


The  term  S(Q  )  can  also  be  separated  into  relative  and  internal 
contributions, 


S(Q°)  =  P°.R°  +  S'(Q'°)   .  (IV. 44) 


With  this  separation  Eq.  (IV. 43)  becomes 


S(Q)  =  P°.R  +  S'(Q'°)  +  [~(V  PR)-dR  +  [~  P'.dQ'   .     (IV. 45) 

Q°~   ~    ~    Q°    ~ 

Since  asymptotically  P_=P„,  we  can  write  the  action  in  the  asymptotic 
region  as 


pQ 
S(Q)  ^^     P°.R  +  S'(Q'°)  +   r  P'.dQ'  (IV. 46) 

Q°"   "     • 

The  only  other  term  to  look,  at  is  \i.      Since  it  contains  the  argument 
of  J  it  will  change  each  time  J  changes  sign.  Asymptotically  this 
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means  it  will  change  by  li  each  time  the  internal  Jacobian  J'  changes 
sign.   The  frequency  of  this  occurrence  depends  on  the  value  of  w. 

To  summarize  the  asymptotic  wave  function,  then,  we  have  the  form 


4:i  -  x°  c. 


r/O 


J' 
a 


i 

e    e 


iii 


P'.dQ'+S 


(Q'°)) 


(IV. 47) 


In  order  to  satisfy  asymptotic  conditions  the  function  in  Eq.  (IV. 47) 

must  be  equivalent  to  the  general  asymptotic  form  in  Eq.  (IV. 30).   Let 

us  first  equate  the  two  functions  yp     L   and  +  at  t  =  t  . 

a,  i 


5  P°.R°   ;  S'(Q'°) 


*i:>(q°,  =  /  C,  e«  '-  •"-  e  « 


(IV. 48) 


^0(Q°)  =  (2iiM)-'^/2  ^^(Q,Oj   ^ 


R  f -^-^ 


(IV. 49) 


Therefore,  setting  \p^   1(0°)  =  +  (Q°)  gives 

a, i  ~       ~ 


X°  Cj.e^  ^'^*^'  ^  =  (2iiM)-^/2  ^^^Q,o^ 


(IV. 50) 


X°Cj,=  (2TiM)-^/2  ^^(Q,Oj  ^-R  S'(^'  > 


(IV. 51) 


•7 

The  factor  C_  is  a  constant  asymptotically.  In  order  for  |C_|  to  be 
the  amplitude  of  the  wave  function  on  electronic  surface  T  it  must  be 
initially  set  equal  to  one.  If  we  put  the  resulting  expression  for  X 
back  into  Eq.  (IV. 47)  we  get 
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*a^r(A)  =  (2nM)  ^'"^   u^(Q'°) 


r/O 


1/2 


X  e 


ml       t     l^i,a<^'L^>l'^^' 


(IV. 52) 


where  the  product  P'.dQ'  has  been  expanded  in  a  sum  over  i  and  dQ' 
"^       ~a  ~  1 


replaced  by  Pi/M  dt. 


,.(-) 


Equation  (IV. 52)  is  the  general  asymptotic  form  of  <|;1  .   It 
satisfies  the  asymptotic  conditions  at  t   (i.e.  at  the  initial  value 
of  R) .  The  conditions  represented  by  Eq.  (IV. 52),  however,  apply  in 
the  entire  asymptotic  region  and  not  just  at  the  initial  value,  of  R 
which  is  chosen  somewhat  arbitrarily  in  that  the  only  prescription  is 
for  it  to  be  outside  the  interaction  region.   Clearly  the  asymptotic 
form  of  ^     does  not  satisfy  initial  conditions  away  from  the  value 
of  R(t  )  due  to  the  terms  involving  J'  and  P' . .   This  is  a  problem 
because  the  transition  amplitude  D^.,  and  hence  the  cross  section, 
should  be  independent  of  the  value  of  R(t  ).   Remember  that  R(t  ) 
refers  to  the  relative  distance  between  the  fragment  X  and  the  atom  A 
after  photodissociation  of  the  molecule.   The  problem  of  satisfying 
asymptotic  conditions  becomes  even  more  complicated  as  the  number  of 
internal  degrees  of  freedom  in  X  increases  (N'  becomes  larger).   The 
problem  is  not  unique  to  this  method,  however,  as  it  is  also  present 
in  purely  classical  calculations  of  transition  probabilities.   An 
approximate  way  of  solving  it  in  such  classical  treatments  is  by 
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taking  suitable  averages  over  the  initial  relative  distance.   We  will 
return  to  this  discussion  in  chapter  VI  where  we  will  show  how  D,.  can 
be  made  to  be  independent  of  R  by  two  different  methods. 

Up  to  this  point  the  treatment  of  the  eikonal  wave  function  has 
been  perfectly  general  in  the  number  of  degrees  of  freedom  allowed  and 
in  the  number  of  excited  states  involved  in  the  molecule  M.  The  only 
constraint  in  this  regard  has  been  that  the  bond  broken  in  the 
photodissociation  event  be  one  between  a  molecular  fragment  X  and  an 
atom  A.  The  feasibility  of  a  more  general  treatment  of  the  problem 
involving  photodissociation  of  a  molecule  into  two  molecular  fragments 
X  and  X'  is  certainly  within  the  scope  of  this  theory  and  should  be 
pursued.   However,  such  a  treatment  is  outside  the  scope  of  this 
dissertation.   Also  outside  the  scope  of  this  dissertation  is  the 
general  solution  to  the  complicated  problem  of  satisfying  asymptotic 
conditions  for  any  number  of  degrees  of  freedom  N.   The  present  work 
has  been  developed  in  terms  general  enough  so  that  such  topics  can  be 
pursued  as  extensions  of  the  basic  formalism  presented  here. 

The  first  step  in  solving  the  problem  of  asymptotic  conditions  in 
general,  as  in  any  investigative  work  is  to  limit  the  discussion  to 
special  model  cases  which  are  more  easily  solved.   In  the  next  section 
we  will  present  various  simple  models  of  photodissociation  so  that  the 
problem  of  meeting  asymptotic  conditions  can  be  solved  in  a  special 
case  and  results  can  be  generated  to  compare  with  previous  results 
from  other  known  theories.   In  this  way  the  accuracy  of  the  theory  can 
be  tested  and  possibilities  for  a  more  general  solution  can  be 
proposed. 


CHAPTER  V 
COLLINEAR  MODELS  OF  POLYATOMIC  PHOTODISSOCIATION 


The  general  photodissociation  event  as  symbolized  by  Eq.  (I.l) 
can  be  greatly  simplified  by  imposing  various  constraints  on  the 
molecule  M  and  the  fragment  X.  The  factors  which  can  be  constrained 
include  the  following: 

i)  the  number  of  internal  degrees  of  freedom  in  X, 

ii)  the  spatial  orientation  of  the  system  throughout  the 
■  process  (i.e.  free  vs.  fixed  orientation), 

iii)  the  number  of  coupled  electronic  states  in  M  (and  X)  involved 
in  the  process. 

Factors  i  and  ii  will  affect  the  number  of  coordinates  in  the  set  ^. 
Factor  iii  will  affect  the  number  of  states  T   to  be  included  and 
consequently  the  size  of  the  <j/  and  *  matrices  in  Eq.  (11.24). 

In  all  of  the  models  presented  in  this  chapter  the  system  will 
have  a  fixed  orientation  in  space  throughout  the  dissociation.   That 
is,  we  will  be  invoking  the  infinite-order  sudden  (lOS)  approximation 
in  dealing  with  the  rotational  degrees  of  freedom  of  the  molecule 
[Pack,  1974].   In  addition  we  will  be  assuming  only  collinear 
dissociation  in  which  all  motion  of  the  molecule  and  the  product 
fragments  occur  along  a  line  defined  by  the  body  fixed  axis  of 
symme  t  ry . 
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Tvo  Electronic  States 

A  simple  model  would  be  one  in  which  there  are  only  two 
electronic  states  r=l,2  (ground  and  excited  states)  involved,  no 
internal  degrees  of  freedom,  and  no  rotation  in  the  fragment  X  or  in 
the  molecule  M  allowed  (i.e.  one  dimensional  motion).  An  example  of 
such  a  model  would  be  the  photodissociation  of  a  single  bond  in  one 
dimension.  The  photodissociation  of  a  single  CI  bond  in  a  polyatomic 
is  such  an  example  and  has  already  been  studied  using  the  Self- 
Consistent  Eikonal  Method  [Swaminathan  and  Micha,  1982].  In  such  a 
case  there  is  one  relative  coordinate  R^n  and  no  internal  coordinates 
Q' .  The  nuclear  wave  function  on  excited  surface  2  simplifies  to 


g/^  \t)   =  (2TtM)  ^^^ 


P° 


P(t) 


^^2   -L{t)       j  ff(t) 

e  ^     e  "     .      (V.l) 


A  slightly  more  complicated  model  is  one  in  which  there  are  again 
two  electronic  surfaces  r=  1  and  2,  but  now  with  one  internal  degree 
of  freedom  in  Q' .   In  this  case  there  are  a  total  of  two  coordinates, 
one  internal  rsQ  and  one  relative  coordinate  R^ri-   This  model 
corresponds  to  the  photodissociation  of  a  linear  triatomic  A — B-C 
where  A  is  the  atom  and  B-C  the  fragment  X.  An  example  of  such  a 
system  is  the  molecule  CH^I  (methyl  iodide).  This  molecule  can  be 
treated  as  a  linear  triatomic  if  one  considers  the  H-  to  move  as  one 
unit  whose  center  of  mass  lies  along  the  C-I  axis  (see  Fig.VII-1). 
This  is  equivalent  to  considering  only  the  umbrella  mode  of  CH-.   It 
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turns  out  that  this  is  a  reasonable  approach  to  CH-I  dissociation 
because  experimentally  it  has  been  shown  that  the  umbrella  mode  of  CH~ 
is  virtually  the  only  one  excited  during  photodissociation  [Shobotake, 
et  al.,  1980]. 

The  set  A  in  this  model  will  consist  of  {t;q'}  where  q'  is  the 
angle  defined  by 


and 


r(t°)  =  A^  sin(q')  3  r° 


Pj.(t°)  =  Ay  0)  cos(q')  a  P° 


(V.2) 


(V.3) 


where  A  is  the  maximum  value  of  r  for  vibrational  level  v,  and  w  is 
V  ' 

the  asymptotic  frequency  of  oscillation  for  u  (r).   The  range  of  q'  is 
0  to  2ii,  which  spans  both  positive  and  negative  values  of  P  (t  )  for 
every  r(t  ).   The  nuclear  wave  function  along  a  trajectory  will  be 


1 
*^"\t;q')  =  (2iiM)^ 


J'°(q') 


J'(t;q') 


2   -2M(t;q')   ^[cT(t;q')-S'(Q'°)] 


Uy(q') 


(V.4) 


where  u  (q')  is  the  value  of  the  harmonic  oscillator  function  at  r(t  ) 
for  vibrational  quantum  number  v.   The  internal  Jacobian  for  this 
model  is  given  by 


J'(t;q')  =  Pj.(t;q')/co 


(V.5) 


The  equation  for  \\i^    '  above  is  exactly  the  same  as  the  general 
form  of  Eq.  (IV. 32)  except  that  the  function  C  is  not  a  matrix,  and 
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the  set  A  is  composed  of  (t;q')'  This  model  will  require  the 
simultaneous  solution  of  8  differential  equations.  There  will  be  four 
equations  for  position  and  momenta,  and  four  to  solve  for  the  terras  Q. 
in  the  Jacobian. 


Three  Electronic  States 

A  more  complicated  model  for  the  photodissociation  of  a 
polyatomic  is  one  which  involves  the  same  number  of  coordinates,  but 
an  additional  electronic  surface  r=3.  Because  we  now  have  two  excited 
potential  energy  surfaces  we  have  to  consider  the  coupling  between 
them.  The  transition  will  be  from  the  ground  state  surface  r=l  to  the 
coupled  surfaces  labeled  by  r=2,  3.   Recall  from  the  transition 
integral  in  Eq.  (IV. 1)  that  mathematically  we  will  have  a  column 
matrix  (2x1)  for  \jr~^   multiplying  a  row  matrix  (1x2)  for  the  dipole 
operator  D  in  the  integrand.  Thus,  the  nuclear  wave  function  on 
surfaces  2  and  3  will  be  the  matrix  of  Eq.  (11.25)  whose  elements  are 
the  expansion  coefficients  for  |Y  (Q)>  with  incoming  wave  boundary 
conditions. 

The  wave  function  for  this  model  will  be 


i<-)(t;q')  =  (2uH)-l/2  e-S'(Q'°)/H 


C(t;q') 


J°(q') 


J(t;q') 


^''^  4u(t;q')   j  a(t;q') 

e  ^       e  "       %(q').  (V.6) 


The  transition  integral  for  the  collinear  model  with  two  degrees  of 
freedom  and  two  excited  electronic  states  can  be  written  as 
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°fi  =  J^*f  ^J^  5-e^  *ilJ|  dR  dr 

=  Jf*2f  521  *  *3f  531I  •  e^  +i  Ml  dR  dr  (V.7) 

where  D^^  and  D^^  are  the  vector  matrix  elements  of  D  and  as  before 
iay  and  ts(^      ,a) . 

The  fact  that  we  have  two  surfaces  brings  up  the  question  of 
which  surface  to  propagate  the  trajectories  on.   Since  the  transition 
can  be  to  either  surface  and  we  want  to  construct  a  wave  function  with 
amplitude  on  both,  then  both  surfaces  must  be  involved  in  the  dynamics 
of  the  problem.   If  we  refer  to  Eqs.  (III. 21)  and  (III. 22)  we  can  see 
how  this  is  possible.   The  equation  for  generating  the  amplitude  part 
of  the  wave  function  (C)  depends  on  the  average  potential  7.   This 
potential  in  turn  depends  on  the  mixing  of  surfaces  2  and  3  by  the 
matrix  C  and  its  adjoint.  Thus  the  matrix  C  is  self-consistently 
coupled  to  the  potential  V.  If  C  C  =  1,  then  we  can  say  that  the 
quantity  |C-[Q(t)]|   represents  the  probability  of  the  system  being  in 
electronic  state  T.      Therefore,  since  the  potential  V  governs  the 
motion  of  a  trajectory,  we  can  say  that  the  electronic  transitions 
between  states  2  and  3  are  self-consistently  coupled  to  the  nuclear 
motion  and  that  ^\~     is  being  propagated  on  the  average  surface  V. 

This  model  will  require  the  simultaneous  solution  of  16 
differential  equations.  Specifically,  there  will  be  four  equations  to 
solve  for  the  real  and  imaginary  parts  of  C;  eight  equations  to  solve 
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for  the  Q^  in  the  Jacobian,  and  four  equations  for  position  and 
momenta.   In  the  next  section  we  will  consider  this  same  model  but  for 
the  case  when  the  two  excited  electronic  states  are  of  different 
symmetry. 

Symmetry  Aspects  and  Cross  Sections 

Having  more  than  one  excited  electronic  state  introduces  an 
additional  complication  if  the  states  are  of  different  symmetry.   For 
example  in  methyl  iodide  it  has  been  shown  experimentally  that  the 
dominant  transitions  are  from  the  ground  electronic  state  to  two 

excited  electronic  states  with  different  symmetry  types.   The  excited 

3       1 
states  referred  to  are,  in  the  notation  of  Mulliken,  the  Q^  and  Q 

states  which  are  of  symmetry  species  A^  and  E  respectively.   In  order 

for  an  electronic  transition  to  be  symmetry  allowed  the  direct  product 

of  the  species  for  initial  and  final  states  with  the  species  for  the 

transition  operator  must  be  totally  symmetric.   The  ground  state  of 

CH-I  is  of  species  A^ ;  therefore,  the  symmetry  selection  rule  requires 

that  for  matrix  elements  D^^  and  D^^  respectively  of  Eq.  (V.7)  to  be 

nonzero 

k^   X  r(D^)  X  A^  =  A^  (V.8) 


E  X  r(Dp  X  Aj  =  A^  (V.9) 


where  r(D')  is  the  symmetry  species  of  the  electric  dipole  moment  D' 
of  the  molecule.  Equation  (V.8)  corresponding  to  a  transition  to  the 
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3 
Q«  state  is  nonzero  only  for  r(D/.)  a  A^  .  According  to  the  character 

table  for  the  point  group  C-   the  species  A^  corresponds  to  D' .   The 

terra  D'  refers  to  the  transition  dipole  operator  along  the  body  fixed 

z'-axis  which  coincides  with  the  symmetry  axis  of  the  molecule.   In 

this  type  of  transition  then,  the  electric  field  must  have  a  component 

oriented  parallel  to  the  body  fixed  z'-axis.   We  will  refer  to  this  as 

the  parallel  transition.   Similarly  Eq.  (V.9)  corresponding  to  a 

transition  to  the  Q  state  requires  that  r(D')  =  E.   According  to  the 

character  table  the  x  and  y  components  of  the  transition  dipole  moment 

(D'  and  D' )  have  E  type  symmetry.   Thus  for  Eq.  (V.8)  to  be  nonzero 
~x     ~y 

the  electric  field  vector  must  have  a  component  perpendicular  to  the 
body  fixed  axis  and  this  is  designated  as  the  perpendicular 
transition. 

Because  of  the  product  D.e  of  Eq.  (V.7)  the  relative  magnitudes 
of  these  transitions  obviously  depends  on  the  orientation  of  the 
molecule  with  respect  to  the  electric  field  vector.   The  factor  D  is  a 
column  matrix  of  vectors  which,  for  the  present  model  of  methyl 
iodide,  consists  of 


D  3 


^-21 


K^31) 


(V.IO) 


D'  =  e  ,D'  +  e  ,D'  +  e  ,D'  (V.12) 

~M   ~x'  X   ~y'  y   ~z'  z  ^    ' 


where  e  ,  (ff'=x',  y' ,  z')  are  the  body  fixed  unit  vectors. 
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In  order  to  carry  out  the  dot  product  in  Eq.  (V.7)  the  elements  of 
Eq.  (V.IO)  will  be  defined  in  terms  of  unit  vectors  as 


521  ^  ^21  ^x' 


(V.13) 


531  '  ^31  ^z' 


where  we  have  arbitrarily  set  the  perpendicular  component  of  the 
transition  dipole  along  the  body  fixed  x'-axis  for  simplicity.  We  are 
allowed  to  do  this  only  because  we  are  dealing  with  a  pseudolinear 
molecule  for  which  the  transition  dipole  is  isotropic  with  respect  to 
rotation  about  the  body  fixed  z'-axis. 

The  direction  of  propagation  of  the  electric  field  in  the 
laboratory  reference  frame  will  define  the  z-axis  of  this  frame. 
Experimentally  this  corresponds  to  the  direction  of  a  laser  beam  which 
is  crossed  at  right  angles  with  a  molecular  beam  of  target  molecules. 
To  obtain  the  experimental  results  that  we  will  be  comparing  with  in 
chapter  VIII  [Shobotake  et  al.,1980;  VanVeen  et  al.,  1984]  the 
researchers  have  used  laser  light  linearly  polarized  in  the  plane  of 
the  crossed  beams.   We  will  define  this  to  be  the  space  fixed  y-z 
plane  as  in  Fig.  (V-1).   With  these  definitions  we  can  use  the  Euler 
angles  to  express  the  dot  product  as 


D.e  = 
~  -ff 


-21  ~a 
k~31  -a) 


21 


'dj-,    cosy  cosP  cosa  -  d„^  siny  sina 
do,  sinp  sina 


(V.14) 
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CH.I  beam 
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>  Z 


Detector 


Fig.  (V-1)  Diagram  of  collision  angles  in  the  spaced  fixed  frame. 
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where  e  is  a  unit  vector  along  the  spaced  fixed  y-axis.   The  angles 

P  and  a  are  the  polar  and  azimuthal  Euler  angles  respectively  of  the 

body  fixed  system  relative  to  the  space  fixed  system.  The  angle  y   is 

the  Euler  angle  for  rotation  of  the  body  fixed  frame  about  the  body 

fixed  z'-axis. 

The  transition  amplitude  for  the  collinear  model  with  two  degrees 

of  freedom  has  been  given  by  Eq  (V.7).   Since  the  dot  product  involves 

a  sum  over  two  states  of  different  symmetry  the  relative  contributions 

from  r'=2  and  r'=3  are,  as  we  have  shown,  dependent  on  the  orientation 

of  the  molecule  with  respect  to  the  direction  of  the  electric  field. 

If  we  invoke  the  lOS  approximation  we  can  allow  the  molecule  to  have 

different  orientations  in  space  and  retain  the  same  form  for  the 

transition  amplitude.  The  value  of  D,.  will,  however,  be 

parametrically  dependent  on  the  Euler  angles  a,  P,  and  y   of  the 

molecule-fixed  coordinate  system  with  respect  to  the  space-fixed 

system.  We  will  denote  this  angle-dependent  integral  in  the  lOS 

(s) 
approximation  as  D^  '  (P,a, y)  where  the  initial  and  final  states 

have  been  specified  in  detail  as  in  Eq.  (11.11).   Note  that  the 

initial  and  final  quantum  states  are  specified  by  ijh(1,v  ,v  )  and 

X  y 

aa(r,v).   By  substituting  Eq.  (V.14)  into  Eq.  (V.7)  we  obtain 


°ap  .,y  (P'°''Y)  =  j['<'2f(R''^»9>'<')  d22(R,r)  (cosycospsina  +  sinycosa) 

^L  1 

+  r|;2£(R,r,P,a)  d^^CRjr)  sing  sina]  <i»^(R,r,  3,  a)  |j  |dRdr  .  (V.15) 
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At  this  point  it  is  necessary  to  define  angles  9  and  <j>  which 
define  the  direction  of  the  detector  in  the  space-fixed  frame.   The 
angles  9  and  <^   are  also  the  polar  and  azimuthal  angles  respectively  of 
the  vector  £,. .  These  angles  are  distinct  from  P  and  a  which  define 
the  orientation  of  the  molecule  at  the  time  of  absorption  of  the 
photon  while  9  and  <i>   define  the  direction  of  detection  for  the  emitted 
fragment.   It  is  important  to  be  very  clear  with  these  distinctions  in 
order  to  properly  interpret  the  quantities  being  calculated.   Recall 
that  in  a  photodissociation  experiment  it  is  the  direction  of  the 
incident  radiation  which  defines  the  space  fixed  z-axis  while  the 
direction  of  detection  is  the  angle  which  the  differential  cross 
section  is  dependent  on. 

In  order  to  account  for  all  possible  orientations  of  the  molecule 

upon  absorption  of  a  photon,  the  square  of  the  transition  amplitude  in 

Eq.  (V.15)  above  must  be  averaged  in  some  way  over  all  initial  angles. 

The  simplest  way  to  do  this  is  to  integrate  over  all  angles  a,  3  and  y 

and  all  rotational  momenta  p  ,  p^,  and  p  using  the  classical  thermal 

a   p      Y 

partition  function  q  for  rotation  and  a  Boltzman  distribution 
function  of  the  classical  rotational  energy  function  H  .   These 
functions  have  the  form  [McQuarrie,  1976] 


q^  =  8ii^Ik_T  (V.16) 

r        o 


Hj.  =  (Pp  +  p^/sin^p)/(2I)  (V.17) 


where  I  is  the  moment  of  inertia  and  k_  the  Boltzman  constant.  The 

o 

momentum  p  does  not  appear  because  we  are  dealing  with  a  linear 
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molecule.   The  appropriate  integral  is 


l°fil?h  =kH'^^^Pe^Pa^/"''^''''V^|,,u^^'°"^>l'   (^-^^^ 

where  g^.  indicates  that  there  is  a  dependence  on  9  and  <i>.      Since 

(s) 
D^.  (P,a, y)  is  independent  of  p  and  p-  in  the  lOS  approximation  the 

integration  over  the  angular  momentum  is  straight  forward  and  results 

in 


|D..  ij.  =  -^  [dodgdY  sin3  |D^^^    (p,a,Y)|^    •        (V.19) 
fi'th   Q^2   J  °^fi.y 

The  collinear  models  introduced  in  this  chapter  describe  one 
dimensional  motion  (motion  confined  to  a  line)  yet  the  perturbing 
influence  (the  electric  field)  is  a  vector  in  three  dimensional  space. 
Note  that  in  the  integral  in  Eq.  (V.19)  we  have  not  yet  invoked  the 
collinear  approximation.   This  is  indicated  by  using  the  vector  g^.  in 
the  subscript  rather  than  the  one  dimensional  scalar  p^..   At  this 
point  we  will  impose  the  restrictions  for  the  collinear  model.   In  a 
collinear  model  we  know  that  the  products  fly  off  along  a  straight 
line  whose  orientation  in  the  spaced  fixed  frame  is  given  by  3  and  a. 
We  can  predict  therefore,  that  the  only  contribution  to  the  cross 
section  amplitude  in  the  direction  of  g^-  will  come  from  an 
orientation  of  (3, a)  which  is  coincident  with  the  g^.  orientation 
given  by  (9,  <t>).  With  these  considerations  we  will  approximate  the 
integral  in  Eq.  (V.19)  by  an  integral  over  the  square  of  the 
transition  amplitude  D^.O, a, y)  from  the  collinear  model.   This 
transition  amplitude  has  two  terms  arising  from  the  perpendicular  and 
parallel  transitions, 
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(2)  (3) 

Dr.Ojtt,  y)    =  Dr.      (cosy  cosP  sina  +   sidy  cosa)   +  Di.    sing  sina 


(V.20) 


where  we  will  define 


°fi'^  =  J'''r'f(R''^)  dj.,;^(R,r)  v(/.(R,r)  |j  |  dR  dr  .       (V.21) 
Using  this  transition  amplitude  we  will  approximate  Eq.(V.19)  by 

l^fi^^'^^'th  =  ^JdedctdY  sing  |D^.O,a,Y)|^  5(cose-cose)  5(^a)  . 
oil 

(V.22) 

Note  that  the  (9,<i>)  angular  dependence  has  been  approximated  by  using 
delta  functions  of  cos9  and  <fr.   The  dependence  of  D^.  on  a  and  g  has 
been  explicitly  determined  in  Eq.  (V.20).   Carrying  out  the 
integration  first  over  y  eliminates  all  cross  terms  giving  the 
equation 

|Dj.(e,<|>)  \l^   =  l^   Jdgda  sing  [  |D^P  |^  (cos^g  sin^a  +  cos^a) 

+  2  Id^J^I^  sin^g  sin^a  ]  5(cose-cosg)  5((^-a)  .   (V.23) 

It  is  the  expression  in  Eq.  (V.23)  that  will  be  integrated  over 
solid  angle  5  to  obtain  the  amplitude  which  will  finally  be  used  to 
calculate  the  integral  cross  section.  The  expression  for  the  integral 
photodissociation  cross  section  within  the  collinear  model  and  in  the 
lOS  approximation,  then,  is 
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2ixM. 
a 


All  that  remains  is  to  carry  out  the  integral  over  9  and  <|).   After 
doing  this  we  arrive  at  an  expression  which  has  the  same  weighting 
factors  for  the  parallel  and  perpendicular  contributions  to  the  cross 
section.   The  final  expression  for  the  integral  cross  section  is 


'H  -  5Kr57r  «<•"'  f  Kfl'-  Kfl')     <^-"' 


This  is  the  expression  we  will  use  in  the  next  chapter  to  calculate 
the  integral  cross  section  for  the  photodissociation  of  methyl  iodide 
using  a  model  with  two  excited  electronic  states. 


CHAPTER  VI 
SATISFYING  ASYMPTOTIC  CONDITIONS 


Statement  of  the  Problem 


At  the  end  of  Chapter  IV  we  began  to  address  the  need  to  impose 
conditions  on  the  asymptotic  form  of  'jy'    •      In  doing  so  we  encountered 
the  problem  that  ;|;   (t;^')  seemed  to  satisfy  asymptotic  conditions 
only  at  the  beginning  of  a  trajectory  (t=t  ).  We  also  hinted  that 
this  problem  would  be  dealt  with  for  a  specific  model  and  a  general 
solution  presented.  In  this  section  we  will  show  how  asymptotic 
conditions  can  be  satisfied  for  a  two  coordinate  linear  model  with  one 
relative  coordinate  R  and  one  internal  coordinate  r.   This  is  the  same 
model  proposed  in  chapter  V.  Two  different  methods  will  be  proposed 
to  accomplish  this.   Each  method  can  in  principle  be  applied  to  any 
two-coordinate  linear  model  regardless  of  the  number  of  coupled 
electronic  surfaces. 

Recall  from  Eq.  (IV. 52)  that  the  asymptotic  form  of  ^         does  not 

seem  to  be  equivalent  to  the  form  of  Eq.  (IV. 30)  except  at  t  .   In 

other  words  ^  t  ^     except  at  the  initial  time.   By  comparing  the 

•^a 

two  functions  and  removing  common  terms  we  arrive  at  the  following 
expression  for  the  internal  part  of  the  asymptotic  wavefunction: 


%(Q°) 


j'° 


1/2   -|y   kl  f    ViV'' 


°         s  u  (Q')  .    (VI. 1) 
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The  root  of  the  problem  is  the  fact  that  u  in  Eq.  (VI. 1)  does  not 
equal  u  ,  where  u  is  the  harmonic  oscillator  vibrational  wavefunction 
for  quantum  number  v.   To  simplify  things  we  will  begin  by  limiting 
the  discussion  to  the  two  dimensional  linear  model.   Within  this  model 
Eq.  (VI. 1)  reduces  to 


Uy(r  ) 


P-(t°.r°) 


P  (t;r") 


1/2  -h  k     l[o    K(^''^°)]dt' 


s  u^[r(t;r")]  . 


(VI. 2) 


Now  that  we  have  identified  the  problem  we  will  proceed  to 

propose  two  methods  to  deal  with  it.  The  first  method  recognizes  that 

u  is  cyclically  equal  to  u  at  time  t=t  and  t=t  +t,  where  x   is  the 

period  of  oscillation  of  u  .   In  the  second  method  we  alter  the  form 
*^  V 

of  ^         to  involve  a  sum  of  two  exponentials  whose  phases  are  equal  in 
magnitude  but  opposite  in  sign.   Using  this  form  we  can  construct  a 
WKB  internal  wavefunction  to  approximate  u  (r). 


Method  1:   Averaging  Over  a  Period 


The  wavefunction  u  (r)  is  a  vibrational  function  for  the  diatomic 

v  ' 

nuclear  vibrational  motion.  There  is  a  frequency  w  and  a  period  x 
associated  with  this  vibrational  motion  where 
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T  =  ^  .  (VI. 3) 


Strictly  speaking  this  period  is  only  present  asymptotically,  but  this 
is  the  only  region  we  are  concerned  with  when  setting  initial 
conditions.  We  know  that  Eq.  (VI. 2)  is  equal  to  u  (r)  at  t  =  t  ; 
therefore,  if  we  can  show  that  Eq.  (VI. 2)  is  also  periodic  over  the 
same  T  then  we  know  that  the  wavefunction  \ji         will  meet  asymptotic 
conditions  at  every  t  =  t  +  nx,  where  n  is  an  integer.   We  will  do 
this  by  showing  that  each  factor  in  u  is  periodic  over  x. 

Si 

The  first  factor  to  consider  is  the  exponential  of  -iy/2.   The 
value  of  y  will  increase  by  2ii  after  each  period  x,  therefore,  the 

exponent  decreases  by  in.  The  other  exponential  factor  contains  the 

2 
integral  of  P'  .  This  integral  increases  by  (2vJi  +  n)  after  a  period 

X,  where  v  is  the  vibrational  quantum  number.   Combining  the  integral 

2 
of  P'   with  vi/2  after  every  period  x  gives  for  the  nth  period. 


=  ni(-ii  +  2vii  +  n)  =  nv2ni   .  (VI. 4) 


From  Eqs.  (VI. 4)  and  (VI. 2)  we  can  conclude  that  the  exponents  in  u 

3. 

add  up  to  an  integral  multiple  of  2lli  after  every  period  x  and,  thus, 
do  not  affect  the  value  of  u  at  times  t=t  +nx  in  the  asymptotic 
region.  The  only  other  factor  involves  the  square  of  the  momenta  and 
will  have  the  same  periodicity  as  the  internal  wavefunction  u  (r) 
because  the  vibrational  momentum  P'(r)  has  the  same  period  x  as  the 
vibrational  motion. 
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We  have  now  looked  at  each  time  dependent  factor  on  the  left  side 
of  Eq.  (VI. 2)  and  have  shown  that  each  returns  to  its  initial  value 
after  every  period  x.      Since  the  function  u  [r(t;r°)]  is  also  periodic 
over  a  time  x,  then,  the  following  equality  holds: 


T 

n 


u^(t  ;r  ) 


P^t°;  r°) 


P  (\ir°) 


l/2-|v.(T^;r°)  i^J   p'(.^;r°)2dt' 


=  u^[r(Tj^;r°)]  (VI. 5) 

where  x  =  t  +nx.  The  equality  in  Eq.  (VI. 5)  shows  that  asymptotic 
conditions  are  met  at  every  x  (n  =  0,1,2...). 

Since  i    satisfies  asymptotic  conditions  periodically,  then  the 
transition  integral  D,.  will  reflect  this  periodicity  and  it  is 
expected  that  D,.  will  be  a  periodic  function  of  the  initial  relative 
coordinate  R(t  )  for  the  bundle  of  trajectories.   Calculations  of  D-. 
with  various  initial  values  of  R  show  this  to  be  true.  This  behavior 
does  not  corresponds  to  physical  reality,  however,  because  the  cross 
section  should  be  totally  independent  of  R(t  ).   Remember  that  R(t  ) 
in  the  two  coordinate  model  is  the  distance  between  A  and  BC  after  the 
photodissociation  event. 

In  order  to  calculate  a  transition  integral  independent  of  R(t  )  s 
R  ,  the  first  method  will  consist  of  simply  averaging  the  square  of 


the  function  D^.(R  )  over  a  distance  in  R  corresponding  to  one  period 
of  oscillat 
defined  by 


9 

of  oscillation  in  r.  This  average  will  be  denoted  Id,.  I    and 

^  I  fi 'ave 
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I-'hIL-nT     T   I"„(R°)|'  (VI.6) 


"  R(T) 


where  Np  is  the  number  of  equally  spaced  values  of  R  in  the  average. 
Note  that  time  t  =0  has  been  arbitrarily  assigned  to  one  end  of  the 
interval.  To  be  consistent  with  a  model  having  incoming  wave  boundary 
conditions  we  will  require  that  R(0)  >  RCt).   The  equation  for  the 
relationship  between  R(0)  and  R(t)  of  Eq.  (VI.6)  is 


R(T)  =  R(0)  +  q-£  (VI. 7) 

"r 


where  P  is  negative. 

From  Eq.  (VI. 7)  we  can  see  that  the  quantity  |D^, |    is 
independent  of  R  .   Thus,  if  we  use  it  in  the  calculation  then  the 
cross  section  will  also  be  independent  of  R  .   For  the  special  case  of 
an  internal  potential  which  is  harmonic  (or  any  symmetric  potential) 
asymptotically,  the  transition  integral  D^.(R  )  is  cyclic  over  only 
half  a  period;  i.e. 


Df.(O)  =  D^.(|)  .  (VI. 8) 

For  such  a  case,  then,  the  average  in  Eq.  (VI.6)  will  be  over  the 
interval  from  R(0)  to  R(t/2). 
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Method  2:   Construction  of  a  WKB  Internal  Wavef unction 


Let  us  look  again  at  Eq.  (VI. 2).   The  idea  behind  the  second 

method  for  meeting  asymptotic  conditions  is  to  construct  ^         in  a  way 

that  would  result  in  y  having  the  form  of  the  WKB  bound  state 

a 

wavefunction. 

The  WKB  wavefunction  for  a  bound  internal  state  has  the  form 
[Messiah,  vol.1], 


"WKB^'^^    = 


J|p'(r)| 


cos 


i-l\/M\ir-l 


-a 


N|p'(r)| 


1  -1 

j|-  S     +     -fj-  S 

e  e 


(VI. 9) 


where  S  =  |J|P  (r)  |dr  -  M  -J 


-a 


=  S 


WKB   A 


k.  =  constant  =  2 


1^ 


amd  M  is  the  reduced  mass  for  motion  along  r. 

Obviously  we  cannot  reproduce  the  sum  of  two  exponentials  in  the 
WKB  wavefunction  unless  we  use  a  sum  of  two  functions  in  \ii       .      Our 
modified  wavefunction,  which  we  will  call  vj^i~  will  have  the  form 
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"Hz'^   =  t  ^  *- 


(VI. 10) 


^  s 


where  ]|/-  =  Xr  e       ,  and   C  =  +>  - 

Asymptotically  we  will  separate  the  action  into  relative  and  internal 
contributions, 


^Ca  =  SRa<^>  -  ^aC^"^^ 


where  the  subscript  "a"  indicates  the  values  of  each  variable  or 
function  in  the  asymptotic  region.   We  can  show  that  this  modified 
wavefunction  satisfies  the  same  differential  equation  as  f       •      Before 
exploring  these  assumptions,  however,  we  will  show  how  ^jT     satisfies 
asymptotic  conditions  for  all  asymptotic  values  of  R. 

Let  us  start  with  the  detailed  asymptotic  form  of  ^jT     so  as  to 
separate  out  that  part  we  can  equate  with  the  WKB  internal 
wavefunction  u™n(J^)'  Asymptotically 


*u 


=  X 


Za 


2   2^a 


IpOR 


K°  e 

+ 


n  a+ 


+  K  e 


i  S 
h  ^a- 


(VI. 11) 


where  we  have  used 


SRa  =PR  • 


(VI. 12) 


Equation  (VI. 11)  must  be  compared  to  \ii    of  Eq.  (IV. 30)  in  order  to 
isolate  u  (r)  which  in  this  case  will  be  equivalent  to  u,n,_(r). 
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Setting  Xj.^ 


(2ii«)  ^^"^   we  obtain 


u^(r)  = 


P  (t°) 


1  1 

2  -  2  ^a 


P^(t) 


K°  e 

+ 


-  S 
M  ^a+ 


I/O 

+  K  e 


\<- 


(IV. 13) 


By  comparing  Eq.  (VI. 5)  with  the  VKB  function  in  Eq.  (VI. 1)  we  have 
the  requirement  that 


lp'(t°)|2  kl  e 


a+ 
h 


+  K  e 


a- 
h 


is      -is 


K^e     +  e 


(VI. 14) 


or 


E  exp 


i(S^^/M  -  M^/2  .  p^) 


=  E  exp 


"JIO 


dr  -  H^ 
4 


/K 


-a 


(VI. 15) 


where  it  is  implied  that 


K°e'''C  =  K° 


K°|p'(t°)|^/2  ^  ^ 


(VI. 16) 
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and   p,  is  a  constant  to  be  determined.  To  prove  the  validity 

of  Eq.  (VI. 15)  and  determine  p,  we  will  first  determine  the  form  of 

S  ,.     To  do  this  we  will  assume  that  S  ,(r)  is  the  asymptotic 

'   o 
internal  action  which  is  calculated  from  a  trajectory  where  P  (r  )  = 


C|P  (rM. 

If  we  define  S,(r  )  as 


r° 


sj.(r°)  =  J  pj.dr  (VI. 17) 


-a 


then  ve  obtain 


S^^(r)  =  J  P^dr  (VI. 18) 


-a 


where  P»  is  the  internal  momentum  whose  value  at  r  is  C|Pr(r  )!•   By 
examining  a  trajectory  whose  initial  internal  coordinate  and  momentum 
are  r  and  |P  (r  ) |  we  can  arrive  at  an  expression  for  S  (r)  in  terms 
of  Syj^(r)  and  A  where 


S^(r)  =  J  |p'(r)|dr  (VI. 19) 

-a 


and  A  is  defined  as  the  number  of  times  that  the  momentum  P  (r)  has 
changed  sign  during  the  trajectory  from  r  to  r.  As  we  follow  this 
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trajectory  we  will  rewrite  the  internal  momentum  in  the  action 
integral  as 


P  (r)  =  ^|P  (r)| 


(VI. 20) 


where  ^  is  the  sign  of  P  .  We  will  also  use  the  relation 


J  |P^(r)  |dr  =  Mii(v  +  1/2) 


(VI. 21) 


where  v  is  the  vibrational  quantum  number.  The  expression  for  S  (r) 
we  arrive  at  is 


S3^(r)  =  I 


S'   (r)  -  ^ 


.4-1 


Ttn  +  7 —  +  A^— 
4      2 


(VI. 22) 


We  can  also  follow  the  trajectory  with  initial  conditions  -r  and 

I  '   o  I  ' 

-|P  (r  )|  to  obtain  S  _(r).   Using  the  same  relations  we  obtain 


S^_(r)  =  I 


S'   (r)  -  ^ 


-l-f 


ilM  +  ^  +  A^  .   (VI. 23) 


The  two  Eqs.  (VI. 22)  and  (VI. 23)  are  close  to  the  form  we  need  to 
construct  the  WKB  internal  wavefunction.  In  order  to  simplify  these 
equations  it  is  important  to  note  the  relationship  of  A  to  ^.   For  Eq. 
(VI. 22)  I,  =   1   for  even  values  of  A  and  -1  for  odd  values  of  A.   This 
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makes  the  sum  (A  +  1/2  -  ^2))i  an  integer  multiple  of  2ii.  Thus,  this 
part  of  S   can  be  dropped  because  the  exponential  of  2iii  is  one.   The 
same  is  true  of  the  sum  (A  +  3/2  -  IJ2)   in  Eq.(VI.23)  because  for  S 
we  find  ^  =  -1  for  even  values  of  A  and  ^  =  1  for  odd  values  of  A. 

At  this  point  it  is  necessary  to  look,  back  at  Eq.  (VI. 15)  and  note 
that  the  asymptotic  form  of  u  is 


M^   =  argfP  (r)]  -  arg[P  (r")]  =  An 


(VI. 24) 


Substituting  Eq.  (VI. 24)  into  the  left  side  of  Eq.  (VI. 15)  and  using 
s'.  from  Eqs.  (VI. 22)  and  (VI. 23)  gives 

3Q 


i^, 


WKB 


4 


=  e 


/M 


iiM 


7  +  P 
4    + 


ij; 


+  e 


WKB 


JiM 
4 


/M 


'WKB 


/M 


-1 


llM 


WKB 


/M 


3ll 


+  P 


(VI. 25) 


where  ^,  is  the  sign  of  P,.   We  are  now  very  close  to  showing  how  the 
two  sides  of  Eq.  (VI. 25)  are  equal.  One  obvious  step  left  is  to 
determine  the  values  of  p  and  p  .   In  order  to  eliminate  the 
exponentials  involving  p^  we  set 


■11/4 


p  =  -311/4  . 


(VI. 26) 
(VI. 27) 


The  final  step  is  to  note  the  relationship  between  I,     and  I,  .     Having 


assumed  that  S   is  calculated  from  the  trajectory  whose  initial 
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momentum  is  P  (r  )  =  |P(r  )|  we  have  prescribed  that  \|;l,~'(t,r  )  will 
be  constructed  from  two  trajectories  whose  position  and  momenta  are 
always  opposite  in  sign  asymptotically;  i.e. 


P_(r)  =  -P^(-r) 


(VI. 28) 


As  a  consequence  of  this  judicious  choice  we  have  that 


^  =  -^. 


(VI. 29) 


Noting  that  cos(2;)  =  cos(-E)  where  Z   is  any  phase  and  using  Eq.(VI.29) 
we  arrive  at  the  equality 


i^. 


nM 


VKB 


/M 


i^ 


+  e 


'WKB 


4 


/M 


1 
tcos 


^VKB^^ 


'      n 
^WKB^"*  "  4 


-1 


^VKB^'^  -  4 


=  e 


+  e 


thus  proving  Eq.  (VI. 15). 

After  so  may  steps  it  is  helpful  now  to  sum  up  the  prescription 
for  this  second  method  of  satisfying  asymptotic  conditions.   We  have 
derived  a  specific  form  for  the  ^,   of  Eq.  (VI. 10) 
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:J;^(t;r  )  = 


2iiM 


2   -  2^ 


I/O 

X  K  e 


(VI. 30) 


where  S,(Q)  =  (T,(t;r  ).  To  satisfy  the  form  of  the  WKB  internal 
wavefunction  asymptotically  it  is  required  that 


K" 


fM~  ^ 

=  2j-^   lp'(t°)|  2 


(VI. 31) 


The  wavefunctions  ^     and  <|(  have  been  defined  as  those  generated  from 
the  trajectories  beginning  at  r  and  -r  with  initial  momenta  |P  (r  ) | 
and  -|P  (r  )|  respectively.  Therefore,  the  function  yj^lT   it',r   )  is 
constructed  from  two  trajectories  by  adding  the  wavefunctions 
generated  along  each  using  the  proper  value  of  p,  from  Eqs.  (VI. 26) 
and  (VI. 27). 

The  wavefunctions  ^     _(t;r  )  of  course  satisfy  the  same 
differential  equation  as  does  ij;   (t;r  ),  because  they  differ  only  by 
the  constant  exp(ip,).   In  this  method  the  transition  integral  of  Eq. 
(IV. 5)  is  independent  of  R  so  there  is  no  need  to  average.  This  is 
due  to  the  way  ftT^   meets  asymptotic  conditions  at  all  (R,r)  in  the 
asymptotic  region.   The  form  of  \j/i   in  this  region  being 
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_1  i, 

,(-),,  o,   ,-,„w,  2  M  P°R 
xp^j.   ^(t;r  )  =  (2tiM)   e 


M 

r 


Jt|P  (t;r°)| 


X  cos 


r 


ij|p'(/)|d.'-f 


-a 


(VI. 32) 


The  major  advantage  in  this  method  over  the  first  one  presented  is 
that  there  is  no  average  over  R  needed  and  thus  not  so  many 
trajectories  are  required.  There  are  two  main  weaknesses  in  this 
approach  however.   The  first  one  is  that  the  internal  wavefunction  is 
being  approximated  by  a  WKB  wavefunction.   This  can  be  a  poor 
approximation  especially  at  low  vibrational  quantum  numbers.   The 
second  weakness  is  the  restriction  that  the  internal  potential  be 
symmetric  asymptotically.   In  most  cases  the  internal  potential  is 
only  approximately  symmetric  asymptotically.  This  particular 
approximation  would  be  especially  poor  at  high  vibrational  quantum 
numbers.   Since  the  WKB  approximation  is  poor  for  small  quantum 
numbers  and  the  symmetric  potential  approximation  poor  at  large 
quantum  numbers  there  will  probably  be  a  certain  range  of  quantum 
numbers  for  which  this  method  gives  optimal  results. 


CHAPTER  VII 
RESULTS:   APPLICATION  TO  METHYL  IODIDE 


In  this  chapter  the  Self  Consistent  Eikonal  Method  is  applied  to 
the  collinear  photodissociation  of  methyl  iodide.   State-to-state 
cross  sections  are  calculated  for  two  different  models.  The  first 
corresponding  to  excitation  from  the  ground  electronic  state  to  a 
single  dissociative  excited  electronic  state.  This  shall  be 
designated  as  the  1-ex  case.  The  second  model  calculation  involves 
excitation  from  the  ground  state  to  two  dissociative  excited 
electronic  states  which  are  coupled  (the  2-ex  case).   Computational 
details  such  as  the  selection  of  initial  conditions  are  also 
presented. 

Coordinates 

Methyl  iodide  is  treated  as  a  linear  triatomic  as  in  the  treatment 
by  Shapiro  and  Bersohn  where  the  umbrella  mode  of  the  CH~  group  is 
modelled  by  an  "equivalent  effective  bond"  resulting  in  a  linear 
triatomic  problem  with  two  vibrational  coordinates.   Essentially,  the 
model  replaces  the  three  hydrogen  atoms  by  an  effective  mass  (of  three 
hydrogens)  located  in  the  plane  of  the  real  hydrogens  and  collinear 
with  the  C-I  bond  axis.  The  coordinates  used  are  the  Jacobi 
coordinates.   In  the  2-coordinate  linear  model  the  Jacobi  coordinates 
correspond  to  the  distance  R  from  the  iodine  atom  to  the  center  of 
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mass  of  CH-  and  the  distance  r  from  the  carbon  atom  to  the  plane  of 

the  three  hydrogens.   In  this  chapter  the  mass  weighted  variables  will 

be  denoted  R,  P,,  and  r,  P  for  relative  and  internal  variables 
K        r 

respectively.  The  distance  R  is  given  in  terms  of  r  and  the 
C-I  distamce  R-,^  as  follows: 


R  =  Rpj  +  0.2011r  .  (VII. 1) 


These  coordinates  are  illustrated  in  Fig.  (VII-1).   Corresponding  to 
each  coordinate  there  is  a  reduced  mass  defined  by 


\   =  (Mj)  (M^jj  )  /  (M^g  j)  (VII. 2) 

Mj.  =  (My  )  (M^)  /  (M^g  )  .  (VII. 3) 


Recall  that  in  chapter  VI  the  problem  of  requiring  ^^~^    to 
satisfy  asymptotic  conditions  at  all  Q  was  addressed  and  solved  by  two 
different  methods.  The  method  of  averaging  over  a  period  was  used  in 
all  the  calculations  for  which  results  will  be  presented  in  this 
chapter.   The  integral  cross  section  expression  for  the  collinear 
model  is  given  by  Eq.  (V.25).   As  discussed  in  chapter  V  the  effect  of 
the  symmetry  of  the  excited  states  on  the  cross  section  is  most 
important  when  considering  two  or  more  excited  states  of  different 
symmetry.   For  the  1-ex  case  the  effect  of  integrating  over  angles  is 
simply  to  multiply  the  cross  section  by  a  constant  1/3.  For  the  2-ex 
case,  however,  the  additional  effect  is  of  have  a  sum  over  the  squares 
of  the  perpendicular  and  parallel  contributions  to  the  transition 
amplitude. 
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Fig.    (VII-1)     Jacobi   coordinates   for  CH-I. 
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When  satisfying  the  asymptotic  conditions  using  method  1  the 
squares  of  the  contributions  are  given  by  the  following  sum 


,(r') 

^fi 


R, 


,(r') 

'fi 


(R°)  r 


R°=  R 


(VII. 4) 


where  the  interval  from  R  to  R,  is  divided  into  N_  equally  spaced 

distances  and  R  >  R.  .   The  value  of  R,  is  chosen  arbitrarily  outside 
a    b  b 

the  interaction  region.  The  length  of  the  interval  (R  -  R,  )  is  such 

that  the  CH»  vibrational  motion  goes  through  half  a  period  as  R  goes 

from  R  to  R,  ,  i.e. 
a     b' 

R  (t  =  t°)  =  R^ 

R  (t°  +  T  /2)  =  R. 


Pr^ 

^b  =  ^a  ^  ~2M^ 


(VII. 5) 


where  x  is  the  period  of  the  umbrella  motion  of  CH-  asymptotically  and 
P„  is  the  incoming  asymptotic  relative  momentum.   Recall  the  variable 
transformation  from  the  set  Q  to  a  set  A  defined  in  Eq.  (III. 2).   In 
the  case  of  the  two  coordinate  linear  model  this  corresponds  to  the 
transformation 


{R,  r}  ->  {t,  w^} 


(VII. 6) 
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where  v^  Is  an  angle  defined  by 


r(t°)  =  A^  sin  w^,  (VII. 7) 


V  is  the  vibrational  quantum  number  for  the  internal  "umbrella"  mode 
of  CH-,  and  A  is  the  value  of  r  at  the  classical  turning  point. 
Because  the  internal  potential  is  asymptotically  harmonic  the  internal 
position  and  momentum  are  given  by 


r(t)  =  A^  sin[w(t  -  t°)  +  w^]  (VII. 8) 


P^(t)  =  Mj.  0)  Aj^  cos  [»(t  -  t°)  +  w^]  (VII. 9) 


P^(t°)  =  Mj.  w  A^  cos  (w^)  (VII.  10) 


where 


(0  = 


(VII. 11) 


Each  trajectory  is  distinguished  by  the  angle  w^ .   A  particular 
value  of  w^  gives  the  initial  internal  position  r(t  )  and  momentum 
p  (t  ).   The  range  of  w^  is  from  0  to  2n.      Thus,  a  grid  of  w-  values 
will  span  all  classically  allowed  values  of  r  along  with  both  the 
positive  and  negative  corresponding  momenta  at  each  r  .  For  example, 
the  value  of  r  has  associated  with  it  both  |P   I  and  -|P   1. 
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This  bundle  of  trajectories  will  all  begin  at  the  same  asymptotic 
value  of  R  with  the  same  initial  time  t  .   As  an  example  of  this  a 


R 


bundle  of  10  trajectories  is  shown  in  Fig.  (VII-2).  There  will  be  N. 
such  bundles  used  in  each  calculation  of  |dL  '^  |  .   For  example  if  the 
range  of  w^  is  divided  up  into  N   equally  spaced  intervals  then  for 
t  =0  there  will  be  N^   trajectories  propagated  from  R  =  R  .   After 

tr  r    I-   o  g 

a  value  for  Id^.  '^(R  )  I  has  been  calculated  another  bundle  of 
'  f 1   ^  a'  ' 

trajectories  is  propagated  from 

R°  =  R(^  )  =  R  and  another  |D|.  (R  )|   calculated.  This 
R 

continues  until  N^  such  bundles  have  been  propagated  and  R,  =  R(t/2) 
is  the  next  R°.   Note  that  [dI.   (Ru)  |   is  not  included  in  the 
average. 

The  initial  relative  momentum  is  given  by  conservation  of  energy 
and  is  negative  indicating  that  the  iodine  atom  is  moving  toward  the 
CH~  fragment.  Remember  that  ;i;    is  being  propagated  in  such  a  way 
that  the  initial  conditions  of  the  trajectories  correspond  to  the 
system  after  the  photodissociation  event.  Thus,  the  system  is 
followed  backward  in  time. 

The  total  energy  of  the  molecule  plus  photon  system  is  conserved 
as  shown  by  the  following 


E^  -  Ej  =  0  (VII. 12) 


p2 

E-  =  -  T^  +  E  +  hv  (VII. 13) 

1     2M     y  ^      ' 

c 


c     R 
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Fig.  (VII-2)  Trajectories  on 
surface. 


the  excited  electronic  potential  energy 
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where  E.  is  the  energy  of  the  CH-I  molecule  plus  photon  system 
initially  (in  the  ground  state),  E^  is  the  energy  of  the  system  after 
having  absorbed  the  photon,  E  is  the  energy  of  the  CH-I  molecule 
initially  for  quantum  state  u=(l,v  ,v  ),  E  is  the  combined  energy  of 
the  fragments  in  the  final  quantum  state  ot=(r,v),  M  is  the  mass  of 
CH-I  and  p^:.sP°       The  numbers  T   and  v  label  the  I  electronic  state  and 

J        1 1   K 

the  CH-  vibrational  state  respectively.  As  in  chapter  II  we  assume 
that  the  momentum  due  to  the  photon  is  much  smaller  than  the  center  of 
mass  momentum  and  so 


P  =  P  +  Mko.  -  P  .  (VII. 15) 

c    c     Ph    c 


Equation  (VII. 12)  then  reduces  to 

P°2 
E  +hv-;^-E  =0  (VII.  16) 

which  gives  the  following  equation  for  P„  : 


P°  =  [2Mj^(E^  +  hv  -  E^)]^^2   .  (VII.  17) 


For  a  given  initial  and  final  quantum  state  E  and  E  are  fixed  and 
the  photon  energy  hv  can  be  varied  from  calculation  to  calculation, 
giving  a  curve  of  cross  section  vs.  photon  energy. 
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One  Excited  Surface 

The  only  states  considered  in  the  first  calculation  (1-ex)  are  the 

3  3 

ground  state  and  Q^  excited  state.  The  Q^  state  leads  to  the 

*  2 
fragments  CH-  +  I  (  Pi /o)  while  the  ground  state  leads  to 

2 
CH-+  I(  Po/o)  *s  shown  in  Figure  (VII-3).  The  potential  surfaces 

given  in  analytical  form  were  determined 'empirically  by  Shapiro  and 

Bersohn  [Shapiro  and  Bersohn,  1980].   The  ground  state  surface  is 


^11  ^^CI'  ""^   =  °e  {exp[-0.899(R^,j  -  4.043)]  -1}^ 

+  J  {k  +  (.1100)  exp[-0.4914(R^j  -  4.043)]} 

X  {r  -  0.6197  exp[-0.4914(R^j  -  4.043)]}^  -D^  -E*     (VII. 18) 

2  3 
k  a  the  force  constant  =  0.0363  e  /a^ 

D  =  0.0874  e^/a^, 
e  0 

E*  =  0.0346  e^/a 


and  the  variable  R„^=(R  -  .2011  r)  is  the  carbon  to  iodine  distance. 

All  energies  are  in  Har trees  and  lengths  in  Bohr  units. 

3 
The  Qq  surface  is  given  by 


V23  (R,r)  =  9.618  exp  (-1.40R) 

+  2.604  exp  (-1.20R  +  0.24r)  +  1/2  (0.0362)r^   .       (VII. 19) 


To  have  meaningful  vibrational  quantum  numbers  for  the  ground 
state  the  R,  r  motion  must  be  transformed  to  normal  modes.   This  was 
done  by  Lee  and  Heller  who  arrived  at  the  x  and  y  modes  given  by  the 
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Fig.  (VII-3)  Diabatic  potential  energy  curves  for  CH3I  at  r=0. 
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following  transformation  [Lee  and  Heller,  1982] 


t  \ 

X 


7.830    -0.1762 
10.6183    4.939  ) 


where 


R  =  Rj,j  -4.043  ao 


R^,^  =  r  -  0.6197  ao 


and  R  =  R  +  .2011  R 


CX 


R 


(VII. 20) 


Lee  and  Heller  have  calculated  the  four  lowest  vibrational 
eigenvalues  for  the  ground  electronic  state.   They  are 


(V,, 

V    ) 

— y- 

(0, 

0) 

(1, 

0) 

(2, 

0) 

(0, 

1) 

E  (a.u.) 
.00399346 
.00643926 
.00885801 
.00949585 


Also  given  by  Lee  and  Heller  is  the  vibrational  wavefunction  for 
the  (0,  0)  state  in  terms  of  harmonic  oscillator  wave  functions  in  the 
normal  coordinates, 


Y^^(x,  y)  =  0.9966110  g^^(x,  y) 


'00 


00^ 


-0.0816282  g^QCx,  y)  -  0.0101739  g^^   (x,  y) 


(VII. 21) 
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where  g    (x,  y)  is  the  product  of  harmonic  oscillator  wave  functions 

for  quantum  numbers  v  and  v  .  The  wavefunction  in  Eq.  (VII. 22)  has 
been  normalized  with  respect  to  the  coordinates  x  and  y.   Before  using 
it  in  the  transition  integral  it  must  be  renormalized  with  respect  to 
the  mass  weighted  Jacobi  coordinates  R  and  r  using  the  determinant  of 
the  rotation  matrix  in  Eq.  (VII. 20).  Results  will  be  shown  for 
excitation  from  this  (0,  0)  level  and  compared  to  those  of  Lee  and 
Heller.   Results  for  excitation  from  the  next  three  vibrational  levels 
(1,0),  (2,0),  and  (0,1)  of  the  ground  state  will  also  be  presented  in 
this  chapter.   For  these  higher  vibrational  levels  a  simple  product  of 
harmonic  oscillator  functions  was  used  for  the  ground  state 
vibrational  wave  function. 

The  equations  of  motion  for  the  1-ex  case  are 


r  =  ^^  =  J-  P  (VII. 22) 

3P    "r  '^ 
r 

R  =  If-  =  i-  P^  (VII. 23) 

9V 
P  =  _  M  =  _  —ex  (VII. 24) 

^     3r     3r 

an     9V 
?„  =  -—=-  -^  (VII. 25) 

3R     3R 


where  V   is  given  by  Eq.  (VII. 19).  The  Jacobian  is  numerically 
calculated  from 


' 
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rS  r 

N 

^R 

P 

r 

J 

X  w^ 

J 

3R 
3w^ 

3r 
3w^ 

-h 

3r 

-\ 

f3R  1 
3w 

ij 

=  p 

1 
1 

=  ^R  qJ  -  \  i 


ip-  -  p^ 


_  pl ex     1 


2 
3  V 


3R3r 


3r 


"^qJ 


(VII. 26) 
(VII. 27) 

(VII. 28) 


2  2 

.,        3^V  ,  3^V  , 

d1       ex  -1  ex  ^1 

"        3r  "  3R3?  ^ 


(VII. 29) 


where  the  following  notation  of  chapter  IV  has  been  used: 


QX  = 


3R_ 

3w, 


Qi  = 


3r_ 
3w, 


(VII. 30) 


0   3w, 


3P 
I 

3w, 


(VII. 31) 


It  is  possible  to  arrive  at  an  approximate  expression  for  the  Jacobian 
by  looking  at  the  asymptotic  form  of  the  terms  in  Eq.  (VII. 26).   The 
expression  for  the  asymptotic  Jacobian  is 


J  =  P_  P  /« 
a    R  r 


(VII. 32) 


In  Fig.  (VII-6)  the  function  J  (dotted  line)  is  compares  to  the  exact 
calculation  of  J  for  a  single  trajectory.  From  this  figure  we  can  see 
that  J  coincides  with  J  until  the  interaction  region,  then  begins  to 

3 


deviate  slightly.  It  is  not  until  the  trajectory  begins  to  exit  the 
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interaction  region  that  the  two  differ  significantly,  and  even  then 

the  difference  is  not  great.   If  this  behavior  of  J  is  typical  of  all 

trajectories  then  it  can  be  assumed  that  for  this  potential  surface 

the  approximate  Jacobian  J  is  sufficient  for  the  calculation  of  the 

a 

transition  integral.  This  assumption  is  shown  to  be  correct  according 

to  trial  calculations  using  J  in  place  of  J  in  which  the  difference 

a 

between  the  results  was  less  than  5%.     We  would  not  expect  this 
behavior,  however  for  an  excited  state  potential  energy  surface  for 
which  the  variation  in  one  of  the  position  variables  along  a 
trajectory  is  large  from  one  trajectory  to  the  next.  For  example  it 
can  be  predicted  by  observing  the  upper  surface  in  Fig.  (VII-4)  that 
the  factor  3R/3w^  in  Eq.  (VII. 26)  will  be  very  close  to  zero 
throughout  a  trajectory.   For  all  of  the  results  in  this  chapter, 
however,  we  have  used  the  exact  form  for  J. 

We  will  be  solving  Eqs.  (VII.22)-(VII.25)  and  (VII.27)-(VII.29) 
along  each  trajectory  making  a  total  of  8  simultaneous  linear 
differential  equations  for  this  model  of  2  coordinates  and  one  excited 
surface.   A  modified  version  of  the  subroutine  DE  was  used  to 
integrate  these  equations.   It  is  based  on  a  Adam's  procedure  and 
follows  a  Predictor-Corrector  sequence  with  variable  step  sizes 
[Shampine  and  Gordon,  1975].   Figure  (VII-2)  shows  some  typical 
trajectories  beginning  at  R,  .  For  all  calculations  R,  was  chosen  to 

0 

be  14.0  a.u.  (7.4  A)  which  is  well  outside  the  interaction  region. 

For  the  cross  section  calculation  of  CH^I  photodissociation  in  the 
1-ex  case  each  state-to-state  cross  section  required  the  following  to 
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obtain  a  value  which  was  converged  to  less  than  5% 

^r  =  50 

where  N   is  the  number  of  trajectories  per  bundle. 

Thus  a  total  of  only  250  trajectories  were  needed  for  each  cross 
section  calculation  to  converge  to  within  an  error  less  than  the 


maximum  statistical  error 


'  100.^  „,^3^' 


l>[250" 


Results  for  the  one  excited  surface  cross  section  are  shown  in 
Figs.  (VII-7)  to  (VII-14).   Figure  (VII-7)  shows  partial  cross  section 
results  from  excitation  of  the  ground  (0,0)  state  to  dissociation  with 
internal  vibrational  state  v,  where  v  goes  from  0  to  8.  These  results 
can  be  compared  to  those  obtained  by  Lee  and  Heller  plotted  in  Fig. 
(VII-8).   The  scale  on  both  figures  is  chosen  so  that  the  peak,  height 
of  the  total  cross  section  is  equal  to  one  as  seen  in  Fig.  (VII-9) 
where  the  total  cross  section  curves  are  compared. 

It  is  clear  from  Figs. (VII-7)  and  (VII-8)  that  the  results  from 
the  Self-Consistent  Eikonal  (SCE)  method  are  in  very  good  agreement 
with  the  "exact"  wave  packet  calculations  of  Lee  and  Heller.  Our 
results  give  the  same  trends  for  the  relative  peak  heights,  showing 
v=4  to  have  the  highest  peak,  followed  in  order  by  v=3,5,2,  and  6. 
The  peak  positions  and  shapes  also  appear  to  be  very  close  to  those  of 
Lee  and  Heller.   A  careful  examination  of  the  peak  positions  of  the 
partial  cross  section  curves  shows  that  all  of  the  curves  from  the  SCE 
method  peak  at  slightly  higher  energies  than  the  curves  from  the 
calculations  of  Lee  and  Heller.  The  average  of  these  differences  is 
found  to  be  approximately  200cm~  •   This  is  also  the  difference  in  the 
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peak  positions  for  the  total  cross  section  curves  of  Fig.  (VII-9).   An 
energy  difference  of  200cm~  is  very  small  compared  to  the  energy  of 
the  photon  (~.5%  of  hv) .   However,  it  is  significant  in  relation  to 
the  vibrational  energy  levels  which  are  separated  by  only  about 
630cm~  .   A  comparison  of  the  values  for  FVHM  values  shows  that  the 
partial  cross  section  curves  for  the  SCE  method  are  slightly  narrower 
(by  about  100-300cm~  )  than  for  the  curves  of  Lee  and  Heller.  This  is 
also  true  of  the  total  cross  section  curves  for  which  the  FWHM  values 
differ  by  SOOcra"''".   Figure  (VII-10)  and  (VII-11)  show  experimental 
results  for  partial  cross  sections  at  two  different  photon  energies 
along  with  the  results  from  the  SCE  method  and  from  Lee  and  Heller. 
These  results  show  that  the  observed  CH^  vibrational  distribution 
peaks  at  v=2  for  photon  wavelengths  X=266nm  and  2A8nm.   The  calculated 
results  differ  at  both  wavelengths,  peaking  instead  at  v=3  and  v=5 
respectively.  The  present  results,  however,  agree  very  well  with 
those  of  Lee  and  Heller  at  these  wavelengths.   Note  that  the  relative 
intensities  of  the  two  theoretical  results  are  preserved  for  each 
figure.  The  disagreement  between  the  theoretical  and  experimental 
results  can  be  attributed  primarily  to  the  potential  energy  surface 
[Shapiro,  1986]. 

Results  for  excitation  from  the  ground  (1,0),  (2,0)  and  (0,1) 
states  to  dissociation  with  internal  vibrational  states  from  v=0  to  8 
are  shown  in  Figs.  (VII-12)  to  (VII-14).   These  results  are  also 
divided  by  the  peak  height  of  the  total  (0,0)  curve.   The  cross 
section  vs.  photon  energy  curves  are  very  interesting  with  respect  to 
the  way  they  reflect  the  nodal  structure  of  the  ground  state 
wavefunction.   Each  of  the  partial  cross  section  curves  in  Fig. (VII- 
12)  from  excitation  of  the  first  vibrationally  excited  (1,0)  ground 
state  show  two  separate  peaks.  This  reflects  the  fact  that  the  (1,0) 
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bound  state  wave  function  has  one  node  in  the  normal  coordinate 
composed  primarily  of  the  I-CH~  stretch.   Of  even  greater  interest  is 
the  structure  of  the  total  cross  section  curve  also  plotted  in 
Fig.  (VII-12).  This  curve  has  a  symmetric  shape  with  two  peaks 
separated  by  3950cm~  and  one  mimima  corresponding  to  an  energy  of 
about  38700cm~  .  The  location  of  the  (1,0)  minima  is  the  same  (within 
an  error  tolerance  of  +50cm~  )  as  that  of  the  (0,0)  peak.   It  lies 
exactly  midway  between  the  two  peaks.  The  peak  at  the  higher  energy 
is  slightly  less  intense  than  the  one  at  lower  energy. 

The  cross  section  curves  in  Fig.  (VII-13)  for  excitation  from  the 
(2,0)  vibrational  state  also  exhibit  the  same  nodal  structure  as  the 
electronic  ground  state  vibrational  wavefunction.   Each  curve  now 
shows  three  peaks  and  two  minima.  The  total  cross  section  curve  has 
two  peaks  of  similar  height  and  a  third,  less  intense,  peak  positioned 
between  them.  The  less  intense  peak  is  at  38150cm~  .  This  is  at 
about  300cm   lower  energy  than  midway  between  the  two  more  intense 
peaks  which  are  at  35300cm~  and  41600cm~  .   Of  the  two  intense  peaks 
the  one  at  higher  energy  is  slightly  more  intense. 

The  cross  section  curves  of  Fig.  (VII-14)  are  for  excitation  from 
the  ground  (0,1)  state.   In  this  state  the  part  of  the  vibrational 
wavefunction  corresponding  primarily  to  I-CH-  stretch,  the  x-mode,  has 
no  nodes  while  that  corresponding  primarily  to  the  CH-  umbrella 
motion,  the  y-mode,  has  one  node.   Each  of  the  curves  in  Fig.  (VII-14) 
has  only  one  peak.  The  total  cross  section  curve  has  one  broad 
asymmetric  peak  which  appears  to  be  a  superposition  of  two  peaks 
positioned  only  1300cm~  apart.   If  this  is  the  case  then  by 
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inspection  we  can  place  the  more  intense  peak  to  be  centered  at  about 
38900cm~  +100cm"  and  the  other  at  37600cm~  +100cm"  .  The  midpoint 
between  the  two  peaks  would  then  be  38250cm~  . 

The  only  other  theoretical  results  available  to  compare  with  for 
excitation  from  the  first  three  vibrationally  excited  levels  of  the 
ground  state  are  from  a  coupled  channels  calculation  by  Shapiro 
[Shapiro,  1981].   Unfortunately  in  generating  these  results  there  was 
an  error  in  the  calculation.  This  error  was  later  corrected  but  only 
the  results  for  the  (0,0)  level  of  the  electronic  ground  state  were 
published  [Gray  and  Child,  1984].  The  original  published  results  do, 
however,  show  the  same  trends  as  the  present  SCE  results.  The  peak 
positions  for  the  partial  absorption  curves  are  close  to  those  for  the 
present  results  and  all  of  the  cross  section  curves  have  the  same 
number  of  peaks  as  the  present  results.  The  major  difference  is  in 
the  distribution  of  the  final  vibrational  levels.  The  results  of 
Shapiro  in  general  show  the  maximum  intensity  to  occur  one  or  two 
vibrational  quanta  lower  than  the  present  results. 

There  are  some  distinctive  features  of  these  curves  that  we  have 
observed  here.   Each  set  of  curves  tends  to  exhibit  the  same  nodal 
structure  as  that  part  of  the  ground  state  vibrational  wave  function 
corresponding  to  the  I-CH~  stretching  motion.   For  example  the  (2,0) 
level  has  a  wave  function  with  two  nodes  in  the  x  direction  and  the 
cross  section  curves  for  excitation  from  this  level  all  have  two 
minima.  Another  feature  is  the  trend  that  the  total  cross-section 
curves  are  centered  at  successively  lower  energies  as  we  go  from  the 
(0,0)  to  (0,1)  levels.  All  three  total  cross-section  curves  for 
(1,0),  (2,0),  and  (0,1)  show  slight  asymmetry  in  that  the  peaks  at 
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higher  energies  are  more  intense  that  those  at  lower  energies.  These 
features  will  be  discussed  again  in  the  last  chapter. 

Two  Excited  Surfaces 

The  second  calculation  involves  excitation  from  the  ground  state 

to  two  coupled  dissociative  states  (2-ex  case).   One  of  the  excited 

3 
states  is  the  Q^   state  as  in  the  1-ex  calculation.  The  other  is  the 

1  2 

Q-  state  which  leads  to  the  fragments  CH„  +  I  (  P,,^).  All  three 

surfaces  (F  =  1,  2,  3)  are  empirical  potentials  by  Shapiro  [Shapiro, 

3 
1986].  The  ground  and  Q^  potentials  are  slightly  different  from 

those  used  in  the  previous  calculation.  The  analytical  expression  for 

the  potentials  are 


,2 


+  D  {exp[-0.87094(R_^  -  4.04326) ]-l}' 


+  J  {k  +  0.11958A  expI-0.4914(R^j  -  4.04326)]) 


x{r  -  0.619702  exp[-0.4914(R^j  -  4.04326)]}  ^   (VII. 33) 


„  .-  ^,  _  0.. 71398  expr-0.38597(R)] 
^22^  '  ^   ~        [1  +   exp[1.5(R  -  4.5)]) 


+  0.82978  exp[-1.5  (R  -  .2011  r)]  +  0.048149  exp[-0.5  (R-.2011r)] 
+  2  k  r^  -  0.034642  (VII. 34) 
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V33  (R»  r)   =  16.972  exp[-1.361(R)l 


+  206.68  exp[-2.217(R  -  0.2011r)]  -  37.286exp[-1.7146(R  -  0.2011r)] 


+  I  k  r^  (VII. 35) 


where  k  =  0.036225,  D®  =  0.0874,  and  E*  =  0.034642. 
Because  we  are  working  in  the  diabatic  electronic  representation  of 
Eq.  (11.31)  there  are  off  diagonal  coupling  elements  in  the  V  block  of 
the  electronic  Hamiltonian  matrix  Eq.  (11.33).   This  matrix  is  real 
and  symmetric,  therefore,  V^^  =  V-™-  The  coupling  as  given  by  Shapiro 
is  only  dependent  on  R, 


V23  =  A  exp[-0.5(R  -  4.2)] 


A  =  2.1  X  10"-^   .  (VII. 36) 


The  relative  strengths  of  the  transition  dipole  matrix  elements  are 
given  also  as  a  function  of  R.  The  equations  for  d-.,^  are  for  the 
parallel  transition 

d^^  =1/(1+  exp[2(R  -  9.8)]  )  (VII. 37) 

and  for  the  perpendicular  transition 

d2i   =  0.48304/(  1  +  exp[2(R  -  9.8)]  )   .  (VII. 38) 
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In  the  2-ex  case  the  final  state  of  the  fragments  is  denoted  not 
only  by  the  internal  vibrational  quantum  number  v  but  also  by  the 
electronic  state  T   of  the  iodine  atom.  Since  the  final  state  of  the 
fragments  determines  initial  conditions  for  the  calculation  there  must 
be  some  way  to  indicate  which  potential  energy  surface  a  trajectory  is 
on  initially.  The  nuclear  motion  wavefunction  provides  the  necessary 
quantity  which  gives  the  amplitude  on  one  or  the  other  potential 
surface.   From  Eq.  (V.6)  we  have 


i^-\t;  w^)  =  (2ilM)  ^^^  Aj.(t;w^) 


1/2     i  ,,   , 
-2M(t;w^) 

e 


-nS'(Q'°)   , 


"H 


i/h  (T(t;w.) 


u  (w, )  e 
v^  1' 

where  we  have  set  C  s  A_  indicating  a  particular  choice  for  the 

asymptotic  value  of  C.   The  matrix  form  of  A  is 


A  = 


^A  ^ 

A 
\  J/ 


(VII. 39) 


From  normalization  we  set 


+  *       * 


(VII. 40) 


and  note  that  d(A  A)/dt  =  0.  Therefore  the  quantity  A_,A_, 
(T'  =  2,  3)  is  the  probability  amplitude  of  the  system  being  in 
electronic  state  T' .      One  of  the  initial  conditions,  then,  is  the 
value  of  A„  and  A_  at  t  =  t  .   For  example,  in  the  CH-.  +  I  (  Po/t) 
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channel  we  have  A^(t  )  =1,  A-  (t  )  =0.   As  an  example  of  the  time 

■k  i(  it 

dependence  of  Aj.,Ap,  we  have  plotted  curves  of  AjA-  and  A-A~  vs  time 

in  Fig.  (VII-12)  for  final  states  corresponding  to  X=266nm,  v=2, 

(v  ,v  )  =  (0,0),  and  w^=4ii/3.   In  Fig.  (VII-15a)  the  chosen  final 
X  y  i 

electronic  state  is  r=2  which  requires  that  A2(t  )=1  and  A-(t  )=0. 
Similarly  Fig.  (VII-15b)  corresponding  to  r=3  requires  that  k^it   )=0 
and  A-(t  )=1.  Notice  that  the  time  required  to  reach  the  interaction 
region  is  slightly  longer  for  the  case  r=3.   This  is  to  be  expected 
since  the  asymptotic  potential  energy  for  r=3  is  larger  than  for  r=2. 
This  figure  shows  that  for  a  typical  trajectory  corresponding  to  final 
state  (r,v)  the  amplitude  A_A-  stays  at  the  same  value  as  at  t  until 

the  interaction  region  is  reached.   In  the  interaction  region  the 

* 
value  A_A-  goes  through  a  sharp  decrease,  then  begins  to  oscillate 

with  a  steadily  increasing  frequency  and  decreasing  amplitude 

fluctuation.  The  period  of  oscillation  ranges  from  about  600  to  200 

a.u.  where  1  a.u.  =  2.42  x  10"  s.  These  oscillations  rapidly  tend 

toward  a  constant  value  for  the  amplitude  which  varies  from  one 

trajectory  to  the  next.   Based  on  an  examination  of  numerous 

trajectories  it  has  been  observed  that  the  value  for  A_  A_  usually 

settles  to  some  value  between  0.6  and  0.8.   According  to  the  time 

scale  in  Fig.  (VII-15)  the  photodissociation  process  takes  place 

within  approximately  0.05  ps.  This  is  in  good  agreement  with  a 

calculated  excited-state  lifetime  of  .07ps  by  Dzvonik  et  al.  [Dzvonik 

et  al.,  1974]. 

The  relationship  between  the  average  potential  and  A^,  is 

implicit  in  the  definition  of  V  of  Eq.  (11.22).   For  the  present 

calculation  we  obtain  the  following: 
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7(Q;  k^,    A^)  =  VJ^^k^k^     +  V23  A^A^ 


•k  -k  -k  -k  _1 


(VII. 41) 


: 


It  is  convenient  to  express  the  complex  quantity  A_,as  (X_,+iY_, ) 

^         mi       ^         ^ 

Thus,  the  effective  potential  can  be  written  as 


7  =  A"^  V  A  (A  "^A) 


+  .S  -1 


[V22  (X2  +  Y^)  +  V33  (X^  +  Y^)  +  2V23  (X2X3  +  Y2Y3)] 

2     2     2     2 
(X^^  +  Y^^  +  X2^  +  Y2^) 


(VII. 42) 


The  differential  equation  for  A  is  the  same  as  that  for  C  in 
Eq.  (III. 21), 


M/i  ^  +  (W)  A  =  0 


(VII. 43) 


where  W  =  V  -  1  V 


"22   "23 

w   w 

I  32  "33J 


(VII. 44) 


giving  the  following  matrix  elements:  W„-=W-„=V^-  and  W__=V_-  V. 
From  Eq.  (VII. 43)  we  get  the  following  for  X_,  and  Y-, : 


^2  "  M  ^^2  ^22  "■  ^3  ^23^  "  3Y2 


(VII. 45) 


^2  "  M  ^^2  ^22  *   ^3  ^23^  "  8X2 


(VII. 46) 
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^3  "  h  ^^2  ^23  "■  ^3  ^33^  "  dY^ 


(VII. 47) 


*3  "  "  M  ^^2  ^23  ^  ^3  ^33^  "  "  9X. 


(VII. 48) 


From  Eqs.  (VII. 45)  to  (VII. 48)  we  can  see  that  X_  and  Y_  are  conjugate 

variables  in  the  Hamiltonian. 

The  differential  equations  for  generating  R,  r,  P„,  P  of  nuclear 

K   r 

motion  and  Q^,  Q^  ,  P^^,  P^  of  the  Jacobian  for  the  2-ex  case  are 


ex 


similar  to  those  for  the  1-ex  case.  The  only  differences  are  that  V 
is  replaced  by  7  and  that  the  time  derivative  of  P.  contains  the 
additional  sum  over  T   due  to  the  dependence  of  V  on  X_,  and  Y_,  as  in 
Eq.  (IV. 19).   The  explicit  form  of  these  time  derivatives  are 


dp; 

dt 


0j> 


C  2-   ^ 
3^V 


X,Y 


l3r8Rj 


^1  ^  3x: 


/  ^ 
97 


X,Y 


l9Rj 


hi 


97 


9Y, 


l9RJ 


^2  ^  9X, 


X,Y 


97 
l9Rj 


X^  +  ^ 
^3   9Y. 


97 


X,Y 


l9Rj 


X,Y 


(VII. 49) 


and 


dp; 

dt 


9  V 


l9R9rJ 


qJ^ 


X,Y 


2-'> 
9  V 


\li^) 


^1  ^  9x: 


X,Y 


97 
l9rj 


X,Y 


3Y. 


97 


l9rj 


yl  ^  L_ 
^2  ^  9X. 


X»Y 


9V 


l9rj 


X^  +  ^ 
^3  *  9Y„ 


X,Y 


9V 
I9rj 


X,Y 


(VII. 50) 
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where  the  factors  X_,  and  Y_,  defined  in  chapter  IV  are  explicitly, 


yl    !i  yl  -  !I2 


1    ^^3  1   ^^3 

^3=3^  ^3  =  9^  •  <^^^-51) 


The  presence  of  these  factors  in  Eqs.  (VI. 49)  and  (VII. 50)  leads  to 
four  additional  differential  equations  to  generate  for  each  X_,  and 
Yj,,  along  the  trajectories.   Equation  (IV. 25)  and  (IV. 26)  give  the 
time  derivative  equations  for  them.   Thus,  for  the  2-ex  case  a  total 
of  16  differential  equations  were  solved  simultaneously  to  generate 
position,  momenta,  and  the  Jacobian  along  each  trajectory. 

To  obtain  results  which  were  converged  to  within  5%  the  2-ex  case 
required  about  the  same  number  of  trajectories  as  the  1-ex  case,  for 
every  state-to-state  cross  section.  This  convergence  was  reached  when 
the  number  of  trajectories  in  each  bundle  was  set  to  50  (N  =50). 
Convergence  in  the  number  of  initial  relative  values  needed  in  the 
average  over  R  was  again  reached  at  N„  =  5.   Thus,  a  total  of  only 
250  trajectories  were  used  for  each  state-to-state  cross  section. 

All  results  for  the  2-ex  case  are  presented  in  Figs.  (VII-16)  to 
(VII-21).   In  Fig.  (VII-16)  all  results  for  the  SCE  method  are  shown 
as  curves  of  cross  section  vs  photon  energy.  The  intensities  for  all 
of  these  curves  have  been  divided  by  the  maximum  intensity  for  the 
total  cross  section.  The  curves  in  Figs.  (VII-16)  and  (VII-17)  are 
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partial  photodissociation  cross  sections  for  final  states  (r,v)=(2,v) 
and  (3,v)  respectively.   In  these  two  figures  the  final  vibrational 
quantum  number  v  for  the  CH~  fragment  ranges  from  zero  to  five.  The 
trend  for  both  r=2  and  r=3  is  for  the  maximum  in  the  vibrational 
distribution  to  shift  from  v=l  to  higher  v  as  the  photon  energy 
increases.  The  most  intense  partial  absorption  for  r=2  is  predicted 
to  occur  at  about  E=42000cm~  and  v=3,  while  for  r=3  the  most  intense 
peak  is  at  an  energy  E=39000cm~  and  for  quantum  number  v=2.   A 
general  feature  of  the  SCE  results  is  for  the  partial  cross  sections 
for  r=2  to  peak  at  energies  which  are  about  2000-2500cm~  higher  than 
those  for  r=3.  The  curves  for  r=2  are  also  much  less  intense  than  for 

r=3. 

Looking  at  the  total  cross  sections  in  Fig.  (VII-18)  which  are 
the  result  of  summing  the  partial  cross  sections  we  see  the  same 
features.  The  total  cross  section  curve  for  r=2  (.)  has  a  peak  at  an 
energy  E„=41900cm~  which  is  2500cm~  higher  than  for  r=3  (+)  which 
peaks  at  E-=39400cm~  .   The  curve  for  the  total  cross  section  has  a 
maximum  at  an  energy  between  E„  and  E-  of  39800cm~  . 

The  features  exhibited  by  these  results  correspond  well  with 
experimentally  known  facts  about  methyl  iodide.  For  instance  it  is 

well  known  experimentally  that  the  transition  to  the  Q  potential 

3 
energy  surface  is  weak  compared  to  that  for  the  Q^  surface  leading  to 

* 
I   [Mulliken,  1940].   The  MCD  absorption  spectrum  for  CH~I  has  been 

estimated  to  be  composed  of  about  25%  intensity  from  Q  and  75Z   from 

3 
Qq  [Gedanken  and  Rowe,  1975].  The  results  from  Fig.  (VII-18) 

indicate  that  for  the  total  photodissociation  cross  section  about  68% 

*  2  2 

shows  up  as  I  (  P-i  /o)  yield  and  32%  as  I(  P-,,^).   The  difference  in 
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the  intensity  distribution  between  the  observed  absorption  spectrum 
and  the  calculated  dissociation  cross  section  is  only  about  9%.   This 
difference  is  actually  to  be  expected  because  the  absorption  spectrum 
and  the  dissociation  cross  section  measure  two  related  but  different 
quantities.  The  relationship  between  the  absorption  and 
photodissociation  cross  sections  will  be  discussed  later. 

In  order  to  help  determine  the  reliability  of  these  results  we 
will  compare  them  with  those  of  Shapiro  [Shapiro,  1986]  from  a  coupled 
channels  type  calculation  using  the  artificial  channel  method 
[Shapiro,  1973].  The  three  graphs  in  Fig.  (VII-19)  are  comparable  to 
those  of  Fig.  (VII-16)  to  (VII-18).   It  is  evident  that  many  of  the 

features  are  the  same.   For  example  about  74%  of  the  total  intensity 

if 
comes  from  the  I  yield.   This  is  more  than  we  predict.   Another 

similar  feature  between  the  present  results  and  those  of  Shapiro  is 

the  trend  for  the  relative  population  of  CH^  vibrational  states  to 

tend  toward  higher  v  as  the  energy  increases. 

The  most  interesting  aspect  of  this  comparison,  however,  are  the 

differences  between  the  two  sets  of  theoretical  results.   The  most 

striking  difference  is  the  peak  position  for  the  total  intensity  of 

2 
I(  Pt/o)  which  is  at  a  much  lower  energy  according  to  Shapiro's 

results  than  for  the  SCE  results.  This  corresponds  to  the  r=3  peak 

for  the  SCE  results.  According  to  Shapiro's  data  the  peak  for  I 

-1  -1  * 

(37700cm  )  is  at  an  energy  550cm   lower  than  the  peak  for  I 

production  (38250cra~  ).   This  is  the  opposite  of  what  we  have  seen 

from  the  SCE  results  where  the  curve  for  I  is  at  a  higher  energy  than 

that  for  I  .  The  peak  position  of  the  r=2  curve,  then,  is  not  at  all 

verified  by  Shapiro's  data  and  represents  a  major  difference  in  the 
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two  theoretical  results  which  cannot  be  explained  simply  by  numerical 
inaccuracies.  The  position  of  the  r=3  peak  is  somewhat  verified  by 
Shapiro's  results,  however.   The  difference  between  Shapiro's  peak,  and 
the  SCE  results  is  not  nearly  as  great  as  for  r=2  but  there  is  still  a 
variance  of  about  1100cm  .  This  is  on  the  order  of  roughly  two 
vibrational  quanta  of  energy. 

The  total  photodissociation  curve  has  a  peak  position  which 
differs  by  much  more  (~1700cm~  )  than  the  r=3  curves  due  to  the  effect 
of  the  curves  for  T=2.     The  total  photodissociation  cross  sections  are 
equivalent  to  total  absorption  cross  sections  because  the  intensity  of 
the  radiation  which  is  absorbed  at  a  particular  photon  energy  must 
translate  directly  into  an  intensity  of  total  fragments  detected.   For 
this  reason  the  curves  for  total  a   can  be  compared  to  the 
experimentally  observed  absorption  curves.   An  absorption  curve  from 

MCD  spectra  [Gedanken  and  Rowe,  1975]  is  able  to  show  the  breakdown  of 

*  2  2 

the  spectrum  in  terms  of  the  transition  to  the  I  (  Pi /o)  ^"<^  I(  ^3/2^ 

surfaces.  This  curve  shows  the  total  absorption  to  peak  at  38500cm~ 

which  is  very  close  to  the  peak  from  Shapiro's  data  at  38100cm~  . 

This  is  to  be  expected  because  Shapiro  fit  his  data  to  these 

absorption  curves  in  order  to  construct  an  empirical  potential  energy 

surface.   The  experimental  partial  absorption  curves  follow  the  trend 

2 
of  the  present  SCE  results  in  that  the  peak  positions  for  I*(  Pi /o) 

2  -12 

and  for  I(  Po/o)  ^^^   separated  by  3300cm   with  the  I(  P-i/n)  peak  at 

higher  energy.  An  interpretation  of  these  different  results  will  be 

given  in  the  next  chapter. 
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We  have  compared  experimental  cross  sections  with  the  present  SCE 
results  and  with  Shapiro's  results  for  X=266nm  in  Fig.  (VII-20)  and 
for  X=248nm  in  Fig.  (VII-21).   Experimental  photodissociation  cross 
sections  for  CH-I  are  only  available  for  photon  wavelengths  of  X=266nm 
and  248nm.   Although  the  curves  of  cross  section  vs  energy  for  the  SCE 
method  differ  substantially  from  Shapiro's  results  the  trends  of  cross 
section  vs  final  v  are  very  similar  for  a  given  photon  wavelength. 
The  intensity  distribution  for  I  yield  vs  v  is  nearly  identical  with 
Shapiro's  at  both  wavelengths.  The  intensity  distribution  for  I  yield 
vs  V  shows  some  discrepancy.   The  SCE  results  for  266nm  and  2A8nm  peak 
at  v=2  and  v=3  respectively  for  the  I  yield  while  Shapiro's  both  peak. 
at  v=2.  The  largest  discrepancy  is  in  the  quantum  yield  at  each 
wavelength.  The  SCE  data  gives  a  yield  for  I  of  8AX   at  266nm  and  67X 
at  248nm.   Shapiro's  data  shows  the  opposite  trend  in  that  the  I 
yield  is  less  (70%)  at  266nm  than  at  248nm  (87%). 

The  comparison  with  experiment  is  much  better  in  the  2-ex 
calculation  than  in  the  1-ex  calculation.  There  are  still  some 
important  discrepancies  however.   At  266nm  the  peak  for  I  yield 
occurs  at  v=2  as  do  both  theoretical  results.  The  peak  for  I  yield  is 
at  v=3  which  is  one  quanta  higher  than  the  theoretical  results.   At 
248nm  the  experimental  data  for  I  and  for  I  yields  both  have  a 
maximum  at  v=2.  The  SCE  results  agree  completely  while  Shapiro's 
results  show  a  maximum  at  v=3  for  I  yield.  According  to  the 
vibrational  distribution,  then,  the  SCE  results  are  slightly  closer  to 
experimental  observation  than  Shapiro's  at  248nm. 
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The  quantum  yields  experimentally  for  both  wavelengths  are  71%. 
As  would  be  expected  the  quantum  yields  from  the  SCE  method  have  the 
opposite  trend  as  those  from  Shapiro's  data.  The  SCE  method  gives 
quantum  yields  of  84%  and  67%  respectively  at  266nm  and  248nm  while 
Shapiro  obtains  70%  and  87%  respectively. 


CHAPTER  VIII 
DISCUSSION  AND  CONCLUSIONS 


In  chapter  II  through  IV  we  have  presented  the  general 
semiclassical  formalism  to  calculate  state-to-state  cross  sections  for 
the  non-adiabatic  photodissociation  of  a  polyatomic  molecule  into  an 
atom  and  a  molecular  fragment.   In  Chapter  V  we  adapted  the  discussion 
from  the  more  general  case  to  the  specific  case  of  a  linear  triatomic 
molecule.  We  discussed  the  wavefunction  and  the  transition  integral 
assuming  collinear  motion  and  invoked  the  lOS  approximation  to  carry 
out  an  average  over  spatial  orientations.   In  Chapter  VI  we  presented 
two  possible  methods  for  satisfying  the  asymptotic  conditions  of  the 
wavefunction  and  developed  them  for  the  model  of  a  collinear 
triatomic.   In  Chapter  VII  we  discussed  in  detail  the  application  of 
the  SCE  method  to  the  collinear  photodissociation  of  methyl  iodide  as 
a  test  of  the  method.  Results  were  presented  along  with  other 
theoretical  results  and  experimental  data  using  first  a  model  with  one 
excited  surface  (1-ex).  To  further  test  the  reliability  of  using  an 
average  potential,  results  were  presented  for  calculations  based  on 
the  same  model  but  with  two  excited  surfaces 

(2-ex).   In  the  first  part  of  this  chapter  we  will  attempt  to  explain 
the  various  features  of  these  results  and  the  deviations  from  other 
theoretical  and  experimental  results  for  the  photodissociation  of 
CH^I. 
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The  1-ex  Case 

The  theoretical  results  presented  in  Fig.  (VII-6)  by  Lee  and 
Heller  are  from  an  exact  time  dependent  wave  packet  calculation  [Lee 
and  Heller,  1982].  These  results  were  later  verified  in  a  revised 
quantum  mechanical  coupled  channels  calculation  by  Shapiro  and  in  two 
semiclassical  calculations  [Gray  and  Child,  1984;  Henriksen,  1985]. 
We  can  conclude,  therefore,  that  the  exact  results  of  Lee  and  Heller 
are  reliable.   In  light  of  this  it  is  most  encouraging  to  find  the  SCE 
results  in  such  close  agreement.   Considering  the  manner  in  which  the 
internal  part  of  the  nuclear  motion  wave  function  for  the  excited 
electronic  state  compares  to  the  expansion  in  harmonic  oscillator 
wavefunctions  used  in  the  exact  treatment,  it  is  very  significant  that 
the  SCE  results  correctly  predicted  the  relative  intensities  for  the 
partial  cross-sections.   Using  harmonic  oscillator  wavefunctions  in 
the  nuclear  motion  wave  function  \j;,  only  to  set  initial  conditions, 
the  SCE  method  correctly  predicts  the  overlap  of  vj;^  with  the 
vibrational  wave  function  t.of  the  ground  electronic  state.   This 
supports  the  assumption  that  a  bound  state  wavefunction  can  be 
adequately  described  by  a  sum  of  exponentials  of  the  action  integral 
along  classical  trajectories  which  span  the  classically  allowed  region 
of  coordinate  space  [Knudson,  Delos,  and  Noid,  1986].  Indeed  many  of 
the  assumptions  made  are  validated  in  this  application. 

Of  particular  importance  is  that  assumption  that  the  most 
important  contribution  to  the  transition  integral  comes  from  the 
classically  allowed  region  of  coordinate  space  and  that  a  finite 
number  of  classical  trajectories  will  adequately  span  this  space. 
This  assumption  is  verified  to  the  extent  that  only  250  trajectories 
were  needed  to  give  converged  results  which  compare  extremely  well 
with  the  exact  results. 
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The  small  differences  which  show  up  in  the  comparison  also  give  us 
information  about  the  effect  of  these  assumptions.  The  small  shift  in 
peak  positions  to  higher  energies  as  well  as  the  smaller  values  for 
FWHM  caul  both  be  explained  by  the  fact  that  the  trajectories,  and  thus 
the  wavefunction,  do  not  extend  into  the  non-classical  region.   The 
contribution  to  the  cross-section  from  the  tail  of  the  quantum 
mechanical  wavef unction  is  not  accounted  for.   According  to  the 
results  this  has  a  greater  effect  on  the  cross-sections  at  lower 
energies  than  it  does  at  high  energies.   The  partial  cross-section 
curves  from  the  SCE  method  fall  off  more  rapidly  on  the  low  energy 
side  than  the  exact  results  do.   If  we  picture  the  overlap  of  ^r   with 
the  wavefunction  t.  we  can  see  that  at  lower  energies  the  quantum 
mechanical  tail  at  small  R  would  overlap  the  vibrational  wavefunction 
in  the  ground  electronic  state  more  than  at  high  energies.   This  is 
also  the  reason  we  see  the  peaks  shifted  to  slightly  higher  energies 
relative  to  the  exact  results.  This  shift  of  200  cm"  practically 
speaking  is  almost  insignificant  if  one  considers  that  the  error  bars 
on  similar  spectroscopic  data  are  estimated  to  be  100  cm"  for  peak 
positions  and  10%  for  curve  widths  (~300cm~  )  [Gedanken  and  Rowe, 
1975]. 

In  comparing  the  theoretical  results  for  the  1-ex  case  with 
experimental  data  at  two  wavelengths  we  find  significant  disagreement. 
Based  on  the  reliability  of  the  theoretical  results  it  can  be 
concluded  that  the  discrepancy  arises  from  either  a  weakness  in  the 
two  coordinates  collinear  model  or  in  the  empirical  potential  of 
Shapiro  and  Bersohn.  As  has  been  mentioned  in  chapter  V  the  only 
vibrational  mode  which  is  significantly  excited  upon  dissociation  of 
CH-I  is  the  umbrella  mode,  which  is  a  symmetric  vibration  parallel  to 


the  axis  of  symmetry  of  the  molecule.  As  a  result  the  collinear  model 
is  expected  to  adequately  describe  the  dynamics  of  dissociation  for 
methyl  iodide.  We  have  also  assumed  that  the  rotational  energy  in  the 
parent  molecule  is  very  small.  Experimentally  the  CH-I  molecules  are 
supersonically  cooled  leaving  rotationally  cold  parent  molecules.   The 
source  of  the  discrepancy  therefore  must  be  the  parametrized  potential 

energy  surfaces.   Shapiro  even  agreed  with  this  conclusion  [Shapiro, 

3 
1986]  and  subsequently  reparametrized  the  Q^  excited  state  surface. 

This  new  surface  is  used  in  the  2-ex  calculations  along  with  a  surface 

for  the  Q  state. 

For  the  1-ex  case  we  also  generated  cross-sections  for  excitation 
from  the  first  three  vibrationally  excited  levels  of  the  ground  state. 
Unfortunately  the  other  theoretical  data  is  in  error  although  the 
general  trends  do  agree.   The  features  of  these  curves  are  very 
reasonable  and  can  be  explained.   As  we  mentioned  in  Chapter  VII  each 
of  the  curves  of  cross-section  vs.  photon  energy  seems  to  exhibit  a 
nodal  pattern  similar  to  that  of  the  ground  state  nuclear  motion 
wavefunction  for  the  normal  mode  x  which  corresponds  primarily  to  the 
I-CH-  stretch.  This  behavior  is  due  to  the  fact  that  different  photon 
energies  correlate  with  different  classical  turning  points  for  the  I- 
CH-j  motion.   Each  photon  energy  in  effect  causes  the  nuclear  motion 
wavefunction  ijj,  to  extend  over  a  different  portion  of  the  ground 
electronic  state  potential  in  the  R-,.,  direction. 

The  trend  for  the  total  cross-section  curves  to  be  centered  at 
successively  lower  energies  is  due  to  the  fact  that  the  energy  gap 
between  the  dissociative  state  and  each  vibrationally  excited  level  of 
the  ground  state  is  successively  smaller.  The  shifts,  however,  are 
generally  about  half  of  what  would  be  expected  from  the  energy 
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difference  between  levels.   It  is  not  clear  why  this  is  the  case.   The 
reason  for  the  asymmetry  of  the  peak  intensities  toward  higher 
energies  is  also  not  clear  but  may  be  due  again  to  the  absence  of  a 
tail  in  the  non-classical  region.  This  asymmetry  is  so  slight, 
however,  that  it  could  just  as  well  be  due  to  numerical  inaccuracy 

The  2-ex  Case 

The  calculations  for  the  2-ex  case  were  different  in  many  respects 

3 
from  the  1-ex  case.   As  was  mentioned  before  the  Q^  surface  was 

reparametrized  by  Shapiro  to  obtain  better  agreement  with  experiment. 

This  calculation  also  involved  another  excited  state  surface  (  Q) 

3 
which,  in  the  diabatic  electronic  representation,  crosses  the  Q^. 

surface.  As  a  result  this  calculation  allowed  for  the  possibility  of 

3  1 

"leakage"  from  the  Q^   surface  to  the  Q  surface  and  vice-verse.   The 

additional  surface  called  for  more  differential  equations  to  be  solved 

(16  compared  with  12  for  the  1-ex  case)  so  there  was  more  CPU  time 

required  per  trajectory.  The  same  number  of  trajectories  were 

required  for  convergence  as  in  the  1-ex  case  yet  about  twice  as  much 

CPU  time  was  required  (~  80s/job  on  the  IBM3080).   A  thermal  average 

over  initial  orientations  and  an  integration  over  detection  angles 

were  also  found  to  be  necessary  for  the  2-ex  case  because  the  excited 

states  are  of  different  symmetry.   For  these  reasons  it  is  difficult 

to  directly  compare  results  for  the  two  cases.   One  similar  feature 

between  the  two  cases  however,  is  the  trend  for  the  maximum  intensity 

to  occur  at  successively  higher  values  of  the  quantum  number  v  as  the 

photon  energy  is  increased. 


As  we  mentioned  in  chapter  VII  for  the  2-ex  case  there  is  a  large 
discrepancy  between  the  results  of  Shapiro  and  the  SCE  results.   The 
best  way  to  begin  to  discuss  this  discrepancy  is  to  compare  both  sets 
of  results  to  the  experimentally  determined  absorption  curves  from  the 

MCD  experiment  [Gedanken  and  Rowe,  1975].  The  absorption  curves 

3      1 
essentially  show  the  relative  amplitude  on  the  Q^  and  Q  electronic 

surfaces  immediately  after  absorption  of  the  photon  but  prior  to 

dissociation  of  the  excited  molecule.  As  the  molecule  is  dissociating 

the  relative  amplitudes  will  be  redistributed  between  the  two  excited 

surfaced  due  to  the  curve  crossing  shown  in  Fig.  (VII-3a).   This  is 

where  the  interesting  dynamics  take  place.   The  redistribution  as  the 

molecule  dissociates  is  obtained  in  the  SCE  method  from  the  behavior 

of  the  amplitudes  A^  and  A^  as  shown  in  Fig.  (VII-12). 

According  to  the  SCE  results  there  is  no  radical  redistribution  of 

amplitudes  as  the  molecule  dissociates.  Evidence  of  this  is  shown  by 

the  peak   for  the  Q  cross-section  curve  which  is  at  about  the  same 

3 
energy  as  the  corresponding  absorption  curve.   The  Q^  cross-section 

curve,  however,  has  its  peak  at  about  1000  cm"  higher  than  the 

absorption  curve  so  there  has  been  some  redistribution  of  amplitudes. 

Further  evidence  of  this  is  given  by  the  difference  in  the  relative 

1      3 
intensity  of  the  Q  and  Q^  curves  for  a  particular  photon  energy. 

For  example  at  a  photon  wavelength  of  266nm  the  experimental 

3      1 
absorption  curves  give  an  intensity  ratio  for  Q^^  to  Q  of  about  20.0 

while  the  SCE  results  give  a  corresponding  cross-section  ratio  of  6.0. 

At  a  wavelength  of  248nm  the  absorption  ratio  is  A.O  while  the  SCE 

cross-section  ratio  is  2.4.   Comparisons  such  as  this  lead  us  to  the 

conclusion  that  there  is  a  substantial  amount  of  "leakage"  from  the 

3  1 

Qj-.  surface  to  the  Q  surface  during  dissociation  caused  by  the  curve 
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crossing.  This  redistribution  is  apparently  greater  at  lower  energies 

and  decreases  at  higher  energies.  At  energies  higher  than  about  41500 

-11  3 

cm  ,  where  the  Q  intensity  becomes  greater  than  the  Q^  intensity, 

this  trend  is  reversed.  The  results  indicate  that  the  "leakage"  at 

-11  3 

energies  above  41500  cm   is  from  the  Q  surface  to  the  Q«  surface. 

We  will  discuss  the  deviations  of  the  SCE  results  from  those  of 

Shapiro  in  terms  of  the  redistribution  of  amplitudes  during 

dissociation.   Shapiro  has  included  partial  absorption  curves  in  the 

published  results.   Because  he  used  the  MCD  results  to  parametrized 

the  potential  energy  surfaces  his  partial  absorption  results  are 

nearly  identical  to  the  experimentally  obtained  curves  of  Gedanken  and 

3 
Rowe.  The  Q^  cross-section  curve  also  peaks  at  about  the  same  energy 

as  the  partial  absorption  curve.  The  cross-section  curve  for  the  Q 

state,  however,  has  a  peak  which  is  at  a  much  lower  energy  than  the 

absorption  peak.  The  difference  is  approximately  4000  cm"  !  The 

ratios  of  cross-section  intensities  at  266nm  and  248nm  are 

respectively  2.3  and  6.7.  This  shows  that  in  Shapiro's  results  there 

3  1 

is  much  more  "leakage"  from  the  Q^.  surface  to  the  Q  surface  at  low 

energies  than  in  the  SCE  results  and  much  less  leakage  at  higher 

energies  than  for  the  SCE  results.   In  fact  already  at  around  an 

energy  of  40000  cm"  the  leakage  trend  has  reversed  and  the  dominant 

1        3 
process  in  the  redistribution  is  for  leakage  from  the  Q  to  the  Q^ 

surface.   Recall  that  this  reversal  also  occurred  in  the  SCE  results 

but  at  an  energy  around  41500  cm"  . 

Both  sets  of  results,  therefore  indicate  a  redistribution  of 


amplitude  from  the  Q^^  surface  to  the  Q  surface  at  low  energies  and 

1  3 

from  the  Q  surface  to  the  Q^  surface  at  high  energies  but  Shapiro's 

results  show  a  much  more  radical  redistribution  than  the  SCE  results. 


Both  sets  of  results  also  show  similar  trends  in  the  vibrational 
distribution  at  each  photon  energy.   At  X  =  266nm  the  distribution  of 
intensity  vs.  vibrational  quantum  numbers  for  both  the  I  and  I* 
channels  from  the  SCE  results  are  identical  to  those  from  Shapiro's 
results  with  respect  to  the  order  of  intensities.   It  is  evident  from 
Fig.  VII-15  that  the  SCE  quantum  yield  is  higher  than  that  of  Shapiro. 
Although  it  does  not  appear  so  by  inspection,  it  turns  out  that 
Shapiro's  quantum  yield  at  266nm  is  almost  identical  to  the 
experimental  value  of  71%.   At  248nm  Fig.  (VII-16)  shows  the  situation 
to  be  reversed  in  so  far  as  the  SCE  quantum  yield  is  close  to  the 
experimental  yield  of  71%  and  Shapiro's  yield  is  much  higher.   From 
this  point  of  view  the  experimental  quantum  yields  at  these 
wavelengths  do  not  greatly  favor  either  set  of  results.   Since  the 
experimental  quantum  yields  are  the  same  at  both  photon  energies  we 
cannot  even  conclude  that  there  is  a  trend  experimentally  for  the 
quantum  yield  to  increase  or  decrease  with  energy.   Other  experimental 
results  are  of  little  help  in  this  regard  as  there  is  not  good 
agreement  between  them  especially  at  248nm  [Riley  and  Wilson,  1972; 
Ozronik  et.  al.,  1974;  Baughcum  and  Leon,  1980;  Hunter  et.  al.,  1978; 
Zhu  et.  al. ,  1987]. 

The  results  from  the  2-ex  case  show  better  agreement  with 

experiment  than  the  1-ex  case.   According  to  Shapiro  this  is  due  more 

3 
to  the  improved  Q^  potential  energy  surface  than  to  the  effect  of  the 

additional  excited  electronic  surface  Q.   The  results  indicate, 

however  that  there  is  still  room  for  improvement  on  the  surfaces. 

Judging  from  Shapiro's  quantum  mechanical  calculations  the  SCE 

method  does  not  seem  to  correctly  model  the  effect  of  the  curve 

crossing  on  the  redistribution  of  amplitudes  as  the  molecule 
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dissociates.   We  can  narrow  the  discrepancy  down  to  the  curve  crossing 
effect  because  the  SCE  results  for  one  excited  surface  are  in  such 
good  agreement  with  the  exact  calculation.  This  failure  to  properly 
account  for  the  curve  crossing  is  very  difficult  to  understand, 
however.   In  calculations  on  a  different  curve  crossing  problem  in 
which  the  collinear  collision  of  H"*"  +  H-  was  studied,  Olson  and  Micha 
found  the  SCE  method  to  give  good  results  compared  with  more  exact 
calculations  [Olson  and  Micha,  1984].   In  this  particular  calculation 
the  transition  integral  was  not  used  and  the  wavef unction  was  not 
explicitly  calculated.   It  was,  therefore,  a  direct  test  for  the 
validity  of  using  an  effective  potential  coupled  to  nuclear  motion  to 
describe  the  dynamics  of  curve  crossing. 

Another  possibility  is  that  Shapiro's  results  are  in  error  with 
respect  to  the  amount  of  redistribution  caused  by  the  curve  crossing. 
The  absorption  step  in  his  calculations  cannot  be  doubted  as  they 
correctly  reproduce  the  MCD  results.   However,  for  dissociation  to  CH- 
+  I  it  is  difficult  to  see  how  there  could  be  such  a  large  shift  in 
the  cross  section  peak  relative  to  the  peak  for  an  optical  transition 
to  this  surface.  As  we  mentioned  previously  this  represents  a  radical 
redistribution  of  amplitudes.   Lacking  other  theoretical  or 
experimental  results  to  validate  Shapiro's  the  possibility  of  error  is 
not  out  of  the  question.   It  should  also  be  noted  that  Shapiro  may 
have  incorporated  some  sort  of  average  over  angular  momentum  states. 
The  formalism  presented  in  his  publication  [Shapiro,  1986]  includes  a 
dependence  on  initial  and  final  angular  momentum  quantum  numbers.   It 
is  not  clear  from  the  article  how  this  was  incorporated  into  the 
results  we  have  been  comparing  to.   If  indeed  they  are  the  result  of 


an  average  over  initial  total  angular  momentum  and  a  sum  over  final 
total  angular  momentum  this  may  have  had  an  affect  on  the 
redistribution  of  amplitudes  which  we  are  not  properly  considering. 

Considerations  on  Angular  Distributions 

In  applying  the  formalism  of  the  SCE  method  in  chapter  VII  we 
used  the  collinear  model  to  calculate  the  cross  section.   In  order  to 
deal  with  the  angular  dependence  of  the  transition  integral  we 
performed  a  simple  thermal  average  over  initial  orientations  of  the 
molecule.  We  also  pointed  out  that  the  cross  sections  were  calculated 
for  the  photodissociation  of  rotationally  cold  molecules.  For  those 
cases  where  the  molecules  initially  have  a  non-negligible  rotational 
energy  it  may  be  important  to  deal  with  the  angular  and  rotational 
states  of  the  system  before  and  after  the  photodissociation  process. 

One  way  of  incorporating  the  angular  momentum  states  in  the 
evaluation  of  the  transition  integral  is  to  begin  with  an  expansion  of 
the  final  state  nuclear  motion  wave  function  in  terms  of  spherical 
harmonics.   For  example  the  wave  function  corresponding  to  an  atom 
diatom  system  could  be  expanded  as 


where  R  and  r  in  the  space  fixed  frame  are  the  vectors  corresponding 
to  the  atom  relative  to  the  center  of  mass  of  the  diatomic,  and  the 
location  of  the  diatomic  atoms  relative  to  each 
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other.  The  solid  angle  co    gives  to  the  direction  of  p^ . 

Efi  ^^'■ 

in  the  space-fixed  frame.  The  quantum  numbers  1  and  m  refer  to  the 
orbital  angular  momentum  and  a  refers  to  the  electronic,  vibrational, 
and  rotational  state  of  the  fragment.   The  next  step  is  to  transform 
the  function  <J;,    to  the  body  fixed  frame  using  rotation  matrix 
elements  D|._.   The  equation  for  this  transformation  is 


'^Ima^^'^-'Pfi)  =  ^   °M,S<°''P''^>  *&'^'X'Pfi>         (^^^^-2) 
J,M,Q 

where  J  and  M  are  the  quantum  numbers  of  the  total  angular  momentum  J 
and  a  is  the  projection  of  J  on  the  body  fixed  axis.  The  angles  a,  P 
and  Y  are  the  Euler  angles  of  the  body  fixed  coordinate  system 
relative  to  the  space  fixed  system  and  X  is  the  angle  between  R  and  r. 
If  this  approach  is  taken  then  the  transition  integral  will  separate 
into  a  sum  of  terms  composed  of  an  integral  over  the  angles  of  (a, P,y) 
and  another  integral  over  the  variables  (R,r,X).   The  first  integrand 
will  contain  a  sum  of  rotation  functions  from  the  final  wavefunction, 
the  angular  dependent  part  of  D.e  ,  and  a  rotational  function  from  the 
initial  state  wave  function.   The  initial  wavefunction  will  be  of  the 
form 


2J.+1  J. 
ti,.(R,r)  =  — i—  Dj^'j^(a,e,Y)  ^^(^,r,X)  (VIII. 3) 

8il      i  i 


where  J.  is  the  initial  total  angular  momentum  quantum  number  and  M. 
and  Q.  are  the  quantum  numbers  for  the  projection  of  J.  on  the  space 
fixed  and  body  fixed  axis  respectively.  The  first  integral  then  will 
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serve  as  a  weighting  factor  for  the  second  integral  which  will  contain 

TMO 

\(;,   ,  i(;.(R,r,X)  and  the  part  of  D.e  dependent  on  (R,r,X)-   The 

TMO 
function  vp,    will  be  calculated  using  the  SCE  method  where  the 

effective  hamiltonian  will  be  dependent  on  J  and  Q. 

An  alternative  to  this  approach  is  to  calculate  the  entire 

function  \j/£(R,r)  including  the  rotational  part  using  the  SCE  method. 

This  will  mean  propagating  the  trajectories  in  terms  of  all  degrees  of 

freedom  (except  the  center  of  mass  motion).  The  transition  integral 

will  be 


°fi  "   H^^'S^  2-tty  '''i(5'E)  ^l   d£  (VIII. 4) 

where  the  vectors  R  and  r  in  the  final  wave  function  appear  in  the 
eikonal  phase  as 


i[S(R°,r°)  +  [p^.dR  +  fp  .dr  ] 
g;.(R,r)  =  X(R,r)  e"  "^  ^      ^    "" 


^fVti't 


(VIII. 5) 


The  integral  will  then  be  expressed  in  terms  of  the  set  A  as 


°fi  =  J*f^^^  5-e^  '"1(6) 


A 


dA  (VIII. 6) 


where  As{t;  9„,  i„,    g^' }  and  0„  and  i>„   are  orientation  angles  of  R  in 
the  spaced  fixed  frame. 

There  is  a  certain  amount  of  freedom  in  defining  c[' .  One 
convenient  set  is  ^'  =  {w^,  Q  ,    <^>  ]   where  9  and  <)>  are  the  angles 
in  w  and  w^  is  the  same  as  in  Eq.  (VII. 7).   The  final  wavefunction 


in  terms  of  A  is 

*f(A)  =  Xf(a')  Cj(A)  e  (VIII. 7) 

where 

X°(a')  =  (2nM)  "^^^  u^(v^)  Yjnj.Cer,*")   •         (VIII. 8) 


V  1   jm.   r  r 
J 


vO   ^0 


The  values  of  9„,  <|>„,  w^,  9  ,  and  <^  will  be  evenly  sampled  in 
K   K   i   r       r 

defining  a  bundle  of  trajectories.   The  magnitude  and  direction  of  the 

initial  incoming  momenta  will  be  set  for  each  trajectory  by 

conservation  of  energy  and  total  angular  momentum.   The  initial  and 

final  total  angular  momentum  must  be  chosen  to  follow  the  selection 

rule  J^  =  J. ,  J.+l. 
f    i'   1- 

Conclusions 

In  this  contribution  we  have  presented  a  semiclassical  formalism 
to  calculate  state-to-state  cross  sections  for  the  photodissociation 
of  polyatomic  molecules  such  as  in  Eq.  (I.l).   This  formalism  employs 
the  calculation  of  a  transition  integral  from  initial  and  final 
wavef unctions.   The  electronic  part  of  the  wavefunction  is  expanded  in 
a  diabatic  electronic  basis  with  coefficients  corresponding  to  nuclear 
motion  wavefunctions.   The  nuclear  motion  wavefunctions  for  the 
electronic  excited  states  are  constructed  in  the  form  of  a  pre- 
exponential  amplitude  matrix  multiplied  by  the  exponential  of  a  common 
phase.  Electronic  transitions  between  coupled  excited  states  during 
dissociation  are  described  using  an  effective  potential  which  is  self 
consistently  coupled  to  nuclear  motion  along  trajectories.   This 


formalism  has  been  shown  to  have  certain  desirable  features  as  stated 
in  the  introduction.  These  features  include  computational  efficiency, 
the  ability  to  deal  with  many  degrees  of  freedom  and  many  coupled 
electronic  states  and  the  ability  to  specify  initial  and  final  quantum 
states. 

In  further  development  of  the  basic  theory  we  have  established  the 
precise  form  of  the  pre-exponential  X  in  terms  of  the  Jacobian  of  a 
transformation  from  the  set  of  coordinates  Q  to  the  set  of  variables 
A.  We  have  developed  the  equations  for  numerically  calculating  this 
Jacobian  exactly  as  a  function  of  time  along  a  trajectory.   We  have 
also  further  developed  the  formalism  by  proposing  two  possible  ways  of 
satisfying  asymptotic  conditions.  As  applications,  we  have  expressed 
the  precise  form  of  the  wavefunction  and  the  transition  integral  for  a 
collinear  model  of  dissociation.   Within  this  collinear  model  we  have 
described  the  angular  dependence  of  the  transition  dipole  for  both 
parallel  and  perpendicular  transitions  and  have  derived  an  expression 
for  the  cross  section  by  performing  a  thermal  average  over  initial 
orientations  of  the  molecule. 

To  test  the  various  aspects  of  this  method  we  performed 
calculations  for  the  dissociation  of  methyl  iodide  within  the 
collinear  model  and  compared  the  results  with  those  from  previous 
theoretical  and  experimental  work.   For  the  calculation  involving  only 
one  excited  surface  we  obtained  excellent  agreement  with  other  well 
established  theoretical  work.  This  calculation  therefore  establishes 
the  validity  of  using  an  eikonal  wavefunction  for  nuclear  motion  on  an 
electronic  surface  having  several  degrees  of  freedom.   It  also 
established  the  validity  of  using  this  eikonal  wavefunction  to 
calculate  transition  amplitudes.  The  small  deviations  from  exact 


results  gave  us  an  indication  of  the  effect  of  ignoring  that  part  of 

the  wavefunction  in  the  non-classical  region,  i.e.  lower  cross 

sections  on  the  low  energy  side.   The  magnitude  of  this  effect, 

however,  was  negligible.   Also  verified  from  the  1-ex  calculation  is 

the  method  of  averaging  over  R  to  satisfy  asymptotic  conditions  as 

well  as  the  method  of  calculating  the  exact  Jacobian  along  a 

trajectory. 

It  should  be  noted  that  in  the  1-ex  case  the  photodissociation 

cross-sections  for  excitation  from  the  excited  vibrational  levels 

(v  ,v  )  of  the  ground  electronic  state  always  required  the  same 
X  y 

bundles  of  trajectories  on  the  upper  surfaces  as  were  used  for  the 

(0,0)  level,  but  with  the  excitation  energies  selected  according  to 

(v  ,v  ).   This  convenient  feature  eliminates  the  need  to  perform  a 
X  y 

completely  separate  calculation  for  each  ground  state  vibrational 
level  and  increases  the  computational  efficiency  of  the  calculation. 
The  results  of  cross  section  vs.  photon  energy  for  excitation  from 
various  excited  vibrational  levels  of  the  ground  state  exhibit  the 
nodal  structure  of  the  ground  state  vibrational  wavefunction  along  the 
dissociation  coordinate. 

As  we  have  discussed  the  results  from  the  2-ex  calculation  are  not 
well  verified  by  the  exact  quantum  mechanical  calculation  of  Shapiro. 
This  seems  to  indicate  that  there  is  a  problem  in  describing  the 
redistribution  of  amplitudes  using  an  average  potential  coupled  to 
nuclear  motion.   When  compared  to  the  experimental  absorption  curves, 
however,  the  SCE  results  are  intuitively  reasonable;  perhaps  even  more 
so  than  Shapiro's  results.  Given  that  there  have  been  no  other 
theoretical  cross  sections  published  to  verify  or  refute  Shapiro's 
results  there  is  a  possibility  that  they  are  in  error.   There  is  a 


need,  then,  for  a  calculation  to  be  performed  by  another  method  using 
these  same  empirical  surfaces  and  couplings  in  order  to  verify  or 
refute  Shapiro's  results  for  the  case  of  two  excited  surfaces. 
Assuming,  however,  that  they  are  correct  it  must  be  concluded  that  the 
aspect  of  the  formalism  which  describes  electronic  transitions  using 
an  effective  potential  must  be  modified. 

The  SCE  results  for  the  2-ex  case  do  show  some  similarities  to 
Shapiro's.  The  trend  of  the  intensities  for  final  vibrational  levels 
is  verified.   Also  verified  by  Shapiro  are  the  overall  magnitudes  of 
the  quantum  yields  for  I*  production  which  hover  around  70%   to  80% 
depending  on  the  photon  energy. 

The  Self-Consistent  Eikonal  method  has  been  shown  to  be  a 
computationally  practical  way  of  calculating  cross  sections  for 
photodissociation  of  polyatomic  molecules.   It  seems  quite  reliable 
for  determining  trends  in  the  partial  cross  sections  and  in  estimating 
the  magnitude  of  quantum  yields.  The  SCE  method  has  been  shown  to  be 
very  accurate  for  photodissociation  cases  with  only  one  excited 
surface. 

The  method  now  should  be  applied  to  other  systems  to  test  other 
features  of  the  theory.   One  extremely  simple  application  would  be 
CD-I  using  the  same  model  potential  energy  surfaces  as  for  CH^I.   This 
would  test  whether  the  theory  correctly  predicts  the  mass  effects  of 
photodissociation.   Another  more  involved  application  would  be  the 
photodissociation  of  CF-I.  The  dynamics  of  this  application  must 
involve  another  internal  degree  of  freedom  besides  the  umbrella 
vibrational  mode.  The  C-F  stretch  vibration  is  very  close  in  energy 
to  the  umbrella  mode,  and  is  probably  significantly  excited  during 
dissociation  [Van  Veen  et  al.,  1985].  It  may  also  be  necessary  to 


include  a  bending  mode  for  the  CF^I  dissociation.   In  that  case  a  full 
three-dimensional  treatment  such  as  we  have  outlined  in  the  last 
section  would  be  necessary.   A  full  three-dimensional  calculation  on 
CF-.I  using  a  wave  function  expansion  has  already  been  performed  by 
Clary  [Clary,  to  be  published].   Once  the  details  of  the  formalism 
have  been  worked  out  for  non-linear  triatomics  it  can  be  extended  to 
other  polyatomics  with  many  degrees  of  freedom.   There  are  a  great 
number  of  polyatomic  systems  for  which  experimental  data  is  available. 
More  experimental  results  are  needed,  however,  in  order  to  establish 
the  reliability  of  theoretical  methods.   It  would  be  particularly 
useful  for  the  present  application  to  have  more  consistent  data  on 
quantum  yields  for  CH-I  photodissociation  as  a  function  of  photon 
energy. 

In  closing  it  is  appropriate  to  summarize  some  of  the  appealing 
aspects  of  the  SCE  method  as  applied  to  photodissociation.   It  is  a 
semiclassical  method  which  employs  classical  trajectories  to  calculate 
nuclear  motion  wavefunctions.  As  such  it  combines  the  desirable 
features  of  classical  calculations  with  those  of  a  quantal  approach. 
From  the  classical  viewpoint  it  provides  insight  into  the  dynamics  of 
the  reaction,  i.e.  motion  of  the  nuclei  during  dissociation  and  the 
time  frame  for  such  motion.  The  use  of  wavefunctions  allow  for 
quantum  mechanical  features  such  as  phase  interference  effects  and 
state-to-state  selection.   Probably  one  of  the  most  appealing  aspects 
of  this  method  over  other  classical  and  semiclassical  methods  is  the 
way  in  which  the  transitions  between  coupled  electronic  surfaces  are 
incorporated.  Unlike  some  other  classical  or  semiclassical  methods 
which  rely  in  part  on  physical  arguments  [Tully  and  Preston,  1971]  the 
expression  for  the  effective  potential  as  well  as  for  the  equation  for 


the  amplitudes  come  directly  from  the  matrix  differential  equation. 
Now  that  this  method  has  been  tested  against  exact  calculations  it  is 
hoped  that  it  will  hold  up  well  against  further  tests  and  become  an 
important  tool  in  predicting  photodissociation  cross  sections  and 
quantum  yields  particularly  in  cases  where  there  are  too  many  degrees 
of  freedom  for  more  exact  methods  to  handle  and  where  several  coupled 
excited  states  are  involved. 


APPENDIX  A 
THE  SHORT  WAVELENGTH  APPROXIMATION 


In  Eq.  (11.58)  we  have  assumed  the  short  wavelength 
approximation.   In  this  approximation  we  require  that 


2M  X 
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7S 


2M 


(A.l) 


or 


«  1 


(A. 2) 


If  we  use  the  Hamilton-Jacobi  form  of  Eq.  (11.51)  we  can  define  a 
reduced  wavelength 


H 


M 


C 


7S 


2M(E-V) 
Thus  we  have  the  requirement  that 


(A. 3) 


.2  ijL 


«   1 


(A. A) 


164 


165 

This  is  the  normal  condition  for  the  short  wavelength  or  classical 
approximation  [Messiah,  I].   The  conditions  of  validity  require  ?t  to 
be  short  compared  with  the  "transverse"  dimensions  of  the  wavefront. 
In  addition  to  this  is  the  requirement  that 

|7X|  «  1   .  (A. 5) 

These  requirements  are  deduced  from  the  optical  analogy  whereby  curves 
of  constant  action  S  form  wavefronts  and  the  rays  drawn  perpendicular 
to  these  fronts  are  parallel  to  VS.  It  is  required  that  the  curvature 
of  the  wavefront  should  be  small  compared  to  1/k. 

If  Eq.  (A. 4)  is  valid  then  we  are  dealing  with  a  case  in  which 
the  wavelength  is  short  and  the  amplitude  X  is  slowly  varying.   It 

turns  out  that  the  amplitude  will  be  roughly  proportional  to  the 

-1/2 
quantity  (P)    .  With  this  in  mind  we  can  make  the  general  statement 

that  the  classical  approximation  calls  for  a  large  but  slowly  varying 

momentum.   Obviously  this  approximation  will  break  down  near  classical 

turning  points  of  one-dimensional  motions  where  the  momentum  goes  to 

zero.   It  will  also  break  down  in  regions  where  the  potential  energy 

surface  changes  rapidly  with  distance  thus  causing  the  momentum  to 

change  rapidly. 


APPENDIX  B 
COMPUTER  PROGRAM  FOR  THE  SCE  METHOD 


Main  Program: 


Call 
ZNLINl 


Call: 

JlNTRP 
)STEP 
FPRIME 


— i. 


Call: 
JCBIAN 
TRAPZ 
PSINV 
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Subroutine  DYNSIM: 


ASIGMA  - 
SIGMAT/NRMAX 


Write 
ASIGMA 


(RETURN 


Do  505 
NR=1,NRMAX 


/\ 


SIGMAT= 
SIGMAT+SIG 


Call 
SIGGEN 


Do  400 
NTRAJ-l, 
NMAN7 


CRTOTl  Z.CRT0T1 +CREAL1 
CIT0T1=CIT0T1+CIMAG1 
CRT0T2=CRT0T2+CREAL2 
CITOT2=.CITOT2+CIMAG2 


Call 
XNLINl 


Call 
DEPHOT 
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The  program  PH0T02  is  the  main  program  vhich  merely  starts  and 
stops  the  program  and  calls  the  driver  DYNSIM. 

The  subroutine  DYNSIM  is  given  in  some  detail  in  the  preceeding 

flowchart.   It  begins  by  reading  from  unit  2  the  namelist  STATES  which 

contains  data  defining  the  final  state  of  the  system.   From  this  data 

all  initial  conditions  are  generated  and  a  DO  loop  over  R  is  entered. 

For  .each  calculation  there  are  NRMAX  bundles  of  trajectories;  each 

with  a  DO  loop  over  NMANY  individual  trajectories.   A  trajectory  is 

propagated  in  the  time  variable  with  DEPHOT,  and  contributes  an  amount 

CREAL  and  CIMAG  to  the  real  and  imaginary  parts  respectively  of  the 

integral.   After  NMANY  such  contributions  have  been  added  the  cross 

section  SIGacrj (R  )  is  calculated  in  SIGGEN.   These  are  averaged  over 

the  R   to  give  ASIGMAscr^.. 
o    *  fi 

The  subroutine  XNLINl  merely  puts  all  initial  values  in  common 
blocks  to  be  used  by  XNLIN2. 

The  subroutine  DEPHOT  is  a  modified  version  of  the  well 
documented  subroutine  DE.   It  is  used  to  integrate  the  differential 
equations.   It  uses  subroutines  INTRP,  STEP,  and  FPRIME. 

The  subroutine  FPRIME  containes  all  of  the  time  derivatives  for 
coordinates,  momenta  and  amplitudes  as  well  as  for  all  Q . ,  P.,  xi  and 
yI  of  the  Jacobian  equations.   It  is  this  subroutine  which  containes 
all  of  the  information  about  the  potential  surfaces. 

The  subroutine  XNLIN2  performs  the  actual  calculation  of  the  real 
and  imaginary  parts  of  the  transition  integral  for  a  single 
trajectory.   It  uses  subroutines  JCBIAN,  TRAPZ,  and  PSINV. 

The  subroutine  SIGGEN  takes  the  real  and  imaginary  parts  of  D^. 
and  computes  the  cross  section  SIG.  The  variable  ASIGMA  is  the  final 
state-to-state  integral  cross  section  for  a  given  photon  energy 
averaged  over  R  . 
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